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Sequence transformations accomplish an acceleration of 
convergence or a summation in the case of divergence by de- 
tecting and utilizing regularities of the elements of the se- 
quence to be transformed. For sufficiently large indices, cer- 
tain asymptotic regularities normally do exist, but the lead- 
ing elements of a sequence may behave quite irregularly. The 
Gaussian hypergeometric series 2-Fi(a, b; c; z) is well suited to 
illuminate problems of that kind. Sequence transformations 
perform quite well for most parameters and arguments. If, 
however, the third parameter c of a nonterminating hyperge- 
ometric series 2-F1 is a negative real number, the terms ini- 
tially grow in magnitude like the terms of a mildly divergent 
series. The use of the leading terms of such a series as in- 
put data leads to unreliable and even completely nonsensical 
results. In contrast, sequence transformations produce good 
results if the leading irregular terms are excluded from the 
transformation process. Similar problems occur also in per- 
turbation expansions. For example, summation results for the 
infinite coupling limit k$ of the sextic anharmonic oscillator 
can be improved considerably by excluding the leading terms 
from the transformation process. Finally, numerous new re- 
currence formulas for the 2^1(0, b; c; z) are derived. 

PACS numbers: 02.30.Gp, 02.30. Mv, 03.65.-w 



s n = s + r n . (1.2) 

The conventional approach of evaluating an infinite series 
consists in adding up so many terms that the remainders 
r n ultimately become negligible. Unfortunately, this is 
not always feasible because of obvious practical limita- 
tions. Moreover, adding up further terms does not work 
in the case of a divergent series since their terms usually 
increase in magnitude with increasing index. 

Alternatively, one could try to determine approxima- 
tions to the remainders r n and to eliminate them from 
the sequence elements s n . At least conceptually, this 
is what a sequence transformation tries to accomplish. 
Thus, the original sequence {snj^o ^ s transformed into 
a new sequence {s^j^Lo wnose elements have the same 
(generalized) limit s but different remainders r' n : 

s' n = s + r' n . (1.3) 

The transformation process was successful if the trans- 
formed remainders r' n have superior numerical properties. 
For example, in the literature on extrapolation methods 
it is said that a sequence transformation accelerates con- 
vergence if the transformed remainders r' n vanish more 
rapidly than the original remainders r n according to 



I. INTRODUCTION 

In mathematics and in the mathematical treatment 
of scientific problems, slowly convergent or divergent se- 
quences and series occur abundantly. Accordingly, many 
techniques for the acceleration of convergence and the 
summation of divergent series have been invented, and 
some of them are even older than calculus (see for in- 
stance pp. 90 - 91 of § or p. 249 of §). 

Sequence transformations are principal tools to over- 
come convergence problems. Let us assume that {s n }^L 
is a slowly convergent or divergent sequence, whose el- 
ements s n may for example be the partial sums of an 
infinite series: 



k=0 



(1.1) 



The basic assumption of all sequence transformations is 
that a sequence element s n can for all indices n > be 
partitioned into a (generalized) limit s and a remainder 
or truncation error r„ according to 
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(1.4) 



Similarly, a divergent sequence {s n }^L 07 whose remain- 
ders r n do not vanish as n — > 00, is transformed into 
convergent sequence {s' n }^L if the transformed remain- 
ders r' n vanish as n — * 00. 

During the last years, considerable progress has been 
reached in this field, as documented by the large number 
of recent monographs p|-|lO[] and review articles jTl"| , |i2| ■ 
Moreover, numerous applications of sequence transfor- 
mations have been reported in the literature. For exam- 
ple, the present author has applied sequence transforma- 
tions successfully in such diverse fields as the evaluation 
of special functions |l2|-[l7j , the evaluation of molecular 
multicenter integrals of exponentially decaying functions 
Jl8|-p2[, the summation of strongly divergent quantum 
mechanical perturbation expansions p^ , |23| -pl[, and the 
extrapolation of crystal orbital and cluster calculations 
for oligomers to their infinite chain limits of stereoregu- 
lar quasi-onedimensional organic polymers |p2|,p3| . 

It should be noted that Pade approximants [p4[ , which 
in applied mathematics and in theoretical physics have 
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become the standard tool to overcome convergence prob- 
lems with power series, can be considered to be a special 
class of sequence transformations since the partial sums 
of a power series are transformed into a doubly indexed 
sequence of rational functions. 

As described above, sequence transformations try at 
least in principle to construct approximations to the ac- 
tual remainders which are then eliminated from the input 
data. This is done by detecting and utilizing regularities 
in the behavior of the elements of the sequence to be 
transformed. For sufficiently large indices n, one can ex- 
pect that certain asymptotic regularities do exist. How- 
ever, sequence transformations are normally used with 
the intention of avoiding the asymptotic domain, i.e., the 
transforms are constructed from the leading elements of 
the input sequence. Unfortunately, sequence elements s n 
with small indices n often behave irregularly. In such a 
case, a straightforward application of a sequence trans- 
formation may be ineffective and even lead to completely 
nonsensical results. Instead, one should analyze the be- 
havior of the input data as a function of the index and ex- 
clude highly irregular sequence elements from the trans- 
formation process if necessary. In this way, one has a 
much better chance of obtaining good and reliable trans- 
formation results. It is the intention of this article to 
describe and classify some of the problems, which can re- 
sult from irregular input data, and to discuss strategies 
to overcome them. 

Sequence transformations are needed most in cases in 
which apart from the numerical values of a few elements 
of a slowly convergent or divergent sequence only very 
little is known. This is a situation which is not uncom- 
mon in scientific applications as for example the sum- 
mation of strongly divergent perturbation expansions as 
they occur quantum mechanics or in quantum field the- 
ory. Thus, it would in principle be desirable to discuss 
problems with irregular input data also via examples of 
numerically determined input data. However, the lack of 
detailed knowledge about the behavior of the elements 
of such a sequence makes it hard to fully understand 
the numerical problems as well as to develop strategies 
to overcome them. Consequently, complications of that 
kind are discussed in this article predominantly via suit- 
able mathematical model problems. 

In Section H, some formal aspect of sequence transfor- 
mations are discussed, in particular the concept of a path 
in the table of a transformation. In this way, different ap- 
proaches for the determination of the approximations to 
the (generalized) limit of the input sequence can be clas- 
sified and formalized. Moreover, the concept of a path 
helps to understand the impact of irregular input data 
on the performance of sequence transformations. 

At least some of the problems mentioned above can be 
illuminated by considering the Gaussian hypergeometric 
function 2-F\(a, b; c; z) which is defined by a power se- 
ries that converges in the interior of the unit circle. This 
function does not only depend on an argument z but also 
on three essentially arbitrary parameters a, 6, and c. As 



discussed in Section III, the convergence of the hyper- 
geometric series can for most values of the parameters 
a, b, and c be accelerated quite effectively by a variety 
of different sequence transformations. Moreover, it is in 
this way frequently possible to associate a finite value to 
a hypergeometric series even if its argument does not lie 
in the interior of the unit circle. 

the situation 



However, as discussed in Section IV 



changes dramatically if the third parameter c of the hy- 
pergeometric series 2-F1 is a negative real number. Then, 
the terms of this series first increase with increasing sum- 
mation index even for \z\ < 1 and produce partial sums 
which look like the elements of a mildly divergent se- 
quence. Only for sufficiently large indices, the terms de- 
crease in magnitude and ultimately produce a convergent 
result. Accordingly, the leading partial sums of such a hy- 
pergeometric series display a highly irregular behavior, 
and they should not be used as input data for a sequence 
transformation. If the leading irregular coefficients are 
skipped and only regular coefficients with higher indices 
are used as input data, then sequence transformations 
are again able to produce good and reliable results. 

Section ^ contains a summary. In Appendix |A|, the 
properties of the sequence transformations, which are 
used in this article, are discussed. Problems with irregu- 
lar input data occur also quite frequently in the mathe- 
matical treatment of scientific problems. In Appendix 
it is shown that the convergence of extensive summation 
calculations for the so-called infinite coupling limit ^3 of 
the sextic anharmonic oscillator can be improved consid- 
erably by excluding the leading irregular coefficients of 
the divergent perturbation series from the transforma- 
tion process. Finally, numerous new recurrence formulas 
for the hypergeometric function 2-Fi(a, 6; c; z) are derived 
in Appendix |^. 



II. ORDER-CONSTANT AND 
INDEX-CONSTANT PATHS 

In this Section, some aspects of sequence transforma- 
tions are discussed which admittedly look very formal. 
Nevertheless, they should not be ignored since they may 
be very consequential in practical applications. 

Obviously, a computational algorithm can only involve 
a finite number of arithmetic operations. Consequently, 
a sequence transformation T can only use finite subsets 
of the original sequence {s„}^ f° r the computation of 
new sequence elements s' m . In addition, these finite sub- 
sets normally consist of consecutive elements. Accord- 
ingly, only subsets of the type {s n , s„+i, . . . , s n+ i} will 
be considered in this article. 

All the commonly used sequence transformations T 
can be represented by infinite sets of doubly indexed 

(n) 

quantities Tf. with fc, n > that can be displayed in 
a two-dimensional array which is called the table of T. 
Here, the convention is used that the superscript n 
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always indicates the minimal index occurring in the finite 

subset of sequence elements used for the computation of 

(n) 

a given '. The subscript k - usually called the order 
of the transformation - is a measure for the complexity 

in) 

of the transformation process which yields T k . 

(n) 

The elements Tf' of the table of T are gauged in 



such a way that T 
sequence element, 



(n) 



corresponds to an untransformed 



T, 



(n) 



(2.1) 



An increasing value of k implies that the complexity of 
the transformation process increases. Moreover, I = l(k) 
also increases. This means that for every k,n > the 
sequence transformation T produces a new transform ac- 
cording to 



T, 



(n) 



MKfe), 



(2.2) 



The exact relationship, which connects k and I, is specific 
for a given sequence transformation T. 

Let us assume that a sequence transformation T should 
be used to speed up the convergence of some sequence 
{ s n} n < = o to its limit s — Sqo. One can try to obtain a 
better approximation to s by proceeding on an in prin- 
ciple unlimited variety of different paths in the table of 
T . Two extreme types of paths - and also those which 
are predominantly used in practical applications - are 
order- constant paths 



rn{n) rp(n+l) rp(n+2) 
1 k y 1 k y 1 k > 



with fixed transformation order k and n 
constant paths 



(2.3) 



oo, and index- 
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(2.4) 



with fixed minimal index n and k — > oo. 

Order-constant and index-constant paths differ signif- 
icantly. It is not even a priori clear that these two types 
of paths lead to the same limit in the case of an arbi- 
trary sequence {s n }" =0 . However, for the sake of sim- 
plicity this potential complication will be ignored here, 
and we shall always tacitly assume that order-constant 
and index-constant paths lead to the same limit. 

In the case of an order-constant path, a fixed number 
of I + 1 sequence elements {s„, s n +i, . . . s n+ i} is used for 
the computation of Tu , and the starting index n of this 
string of fixed length is increased successively until ei- 
ther convergence is achieved or the number of available 
elements of the input sequence is exhausted. 

In the case of an index-constant path, the starting in- 
dex n is kept fixed at a low value (usually n = or 
n = 1) and the transformation order k is increased and 
with it the number of elements contained in the subset 
{s n , s„_|_i, . . . s Jl+ ;( fc )}. Thus, on an index-constant path 
it is always tried to compute from a given set of input 



(n) 

data that element Tj, which has the highest possible 
transformation order k. 

In order to clarify the differences between order- 
constant and index-constant paths, let us consider the 
probably best known sequence transformation, Wynn's 
epsilon algorithm 139]: 
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0, 
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(2.5a) 
(2.5b) 



Wynn |35[ showed that if the input data s„ for the epsilon 
algorithm are the partial sums 



fn{z) = ^2luZ" 
v=0 



(2.6) 



of a (formal) power series for some function f(z), then 

(n) 

the elements e 2j J with evi 
imants to / according to 



(n) 

the elements e 2j J with even subscripts are Pade approx 



>) 

l 2k 



k/k] 



(2.7) 



Here, the notation of the monograph by Baker and 
Graves-Morris Q is used, i.e., a Pade approximant [l/m] 
corresponds to the ratio of two polynomials Pi{z) and 
Qm{z), which are of degrees / and to, respectively, in 

(n) 

z. In contrast, the elements e^fe+i with odd subscripts 
are only auxiliary quantities which diverge if the whole 
process converges. 



It follows from (2.7) that the epsilon algorithm (2.5) 
effec ts the following transformation of the partial sums 
( |2.6|) to Pade approximants: 

{fn(z),.fn+l(z),...Jn+2k(z)} — > [n + k/k}. (2.8) 

Thus, if we use a window consisting of 2k + 1 partial 
sums f n +j(z) with < j < 2k on an order-constant path 
and increase the minimal index n successively, the ep- 
silon algorithm produces the following sequence of Pade 
approximants: 



[n + k/k], [n + k + 1/fc], . . . , [n + k + m/k], 



(2.9) 



Only 2k + 1 partial sums are used for the computation of 
the Pade approximants, although many more are known. 
Obviously, the available information is not exploited op- 
timally on such an order-constant path. 

Moreover, the degree of the numerator polynomial of a 
Pade approximant [n + k + m/k] increases with increasing 
to > 0, whereas the degree of the denominator polyno- 
mial remains fixed. Thus, these Pade look unbalanced. 
Instead, it seems to be much more natural to use diagonal 
Pade approximants, i.e., Pade approximants with numer- 
ator and denominator polynomials of equal degree, or - 
if this is not possible - to use Pade approximants with 
degrees of the numerator and denominator polynomials 
that differ as little as possible. 
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This approach has in principle many theoretical as well 
as practical advantages. For example, Wynn could show 
that if the partial sums fo(z), fi(z), ■ ■ •, f2n{z) of a Stielt- 
jes series are used for the computation of Pade approxi- 
mants, then the diagonal approximant [n/n] provides the 
most accurate approximation to the corresponding Stielt- 
jes function /(z), and if the partial sums fo(z), fi(z), ■ ■ •, 
f2n+i(z) are used for the computation of Pade approxi- 
mants, then either [n+l/n] or [n/n+1] provides the most 
accurate approximation (Theorem 5 of |56|). A detailed 
discussion of Stieltjes series and their special role in the 
theory of Pade approximants can for instance be found in 
Section 5 of the monograph by Baker and Graves-Morris 

& 

Thus, it is apparently an obvious idea to try to use 
either diagonal Pade approximants or their closest neigh- 
bors whenever possible. Let us assume that the partial 
sums fo(z), fi(z), . . ., fm(z) are known. If m is even or 
odd, m = 2/i or m — 2/i + 1 , respectively, the elements of 
the epsilon table with the highest possible transformation 
orders are given by the transformations 

{/o(4/iW, • • • , h„{z)} — > 42 = [/VM] , (2-10) 
{/i(*),/2(*),...,/ 2m +i(*)} — = [m + Vm]- (2.11) 

With the help of the notation [x] for the integral part of 
x, which is the largest integer v satisfying v < x, these 
two relationships can be expressed by a single equation 
(Eq. (4.3-6) of @): 

{ fm—2 

= [m-[m/2]/Im/2]]. (2.12) 

For m = 0, 1, 2, . . ., these transformations correspond to 
the following staircase sequence in the Pade table (Eq. 
(4.3-7) of @): 

[0/0], [1/0], [1/1],... 

...,[v/v],[v + \/v],[v+l/v+ll... . (2.13) 

This staircase sequence exploits the available information 
optimally if the partial sums f m (z) with m > are com- 
puted successively and if after the computation of each 
new partial sum the element of the epsilon table with the 
highest possible even transformation order is computed. 
Moreover, the Pade approximants obtained in this way 
look balanced since the degrees of their numerator and 
denominator polynomials differ as little as possible. 

The example of Wynn's epsilon algorithm strongly in- 
dicates that index-constant paths are at least in prin- 
ciple computationally more efficient than order-constant 
paths since they exploit the available information opti- 
mally. This is in general also true for all other sequence 
transformations considered in this article. 

Another serious disadvantage of order-constant paths 
is that they cannot be used for the summation of di- 
vergent sequences and series since increasing n in the set 



{s n , s„+i, . . . , s n +;} of input data normally only increases 
divergence. 

In view of the examples given above, it is apparently 
an obvious idea to use exclusively index-constant paths, 
and preferably those which start at a very low index n, 
for instance at n = or n = 1. This is certainly a good 
idea if all elements of the input sequence contain roughly 
the same amount of useful information. If, however, the 
leading terms of the sequence to be transformed behave 
irregularly, they cannot contribute useful information, or 
- to make things worse - they contribute wrong informa- 
tion. In such a case it is usually necessary to exclude the 
leading elements of the input sequence from the trans- 
formation process. Thus, it is preferable to use either 
an order-constant path or an index-constant path with a 
sufficiently large starting index n. The use of an order- 
constant path has the additional advantage that the di- 
minishing influence of irregular input data with small in- 
dices n should become obvious from the transformation 
results as n increases. 

Finally, an important theoretical advantage of order- 
constant paths should also be mentioned. Normally, it 
is much easier to perform a theoretical convergence anal- 
ysis of a sequence transformation on an order-constant 
path than on an index-constant path. As the starting in- 
dex n of the string {s n , s„+i, . . . , s„+z} of input data be- 
comes large, asymptotic approximations to the sequence 
elements s n can be used. Often, this greatly simplifies a 
theoretical analysis. In the case of index- const ant paths, 
such a simplification is not possible because not all input 
data have large indices. Accordingly, theoretical conver- 
gence properties of sequence transformations are studied 
almost exclusively on order-constant paths. Notable ex- 
ceptions are two articles by Sidi |37],[58) where the con- 
vergence properties of sequence transformations on both 
order-constant and index- const ant paths are analyzed. 

III. THE GAUSSIAN HYPERGEOMETRIC 
FUNCTION 

The Gaussian hypergeometric function 2-F\(a, b; c; z) is 
one of the most important special functions of mathe- 
matical physics, and its properties are discussed in nu- 
merous books, for example in those by Abramowitz and 
Stegu n pg] , Erdelyi, Ma gnus, Oberhettinger, and Tri- 
comi [|40|| , Magnus, Oberhettinger, and Soni @, Seaborn 
@, Slater Q, Spanier and Oldham @, Temme @, 
and Wang and Guo Q . 

The hypergeometric function 2^1 (a, b; c; z) is defined 
via the corresponding hypergeometric series (p. 37 of [El) 

2 F 1 (a,b;c;z) = ^ ^ z , (3.1) 

m— 

where (a) m = T(a + m) IT{a) is a Pochhammer symbol 
(see for example p. 3 of EJ). The series ( fTl| ) terminates 
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after a finite number of terms if either a or b is a negative 
integer. Otherwise, it converges in the interior of the unit 
circle, i.e., for \z\ < 1, and it diverges for \z\ > 1. On 



the boundary \z\ = 1 of the unit circle, the series (3.1) 
diverges if Re(a + b — c) > 1, it converges absolutely for 
Re(a + b — c) < 0, and it converges conditionally for 
< Re(a + b — c) < 1, the point z = 1 being excluded. 

Thus, a nonterminating hypergeometric series does not 
suffice for the computation of the corresponding hyper- 
geometric function 2 Fx{a, b; c; z), which is in general a 
multivalued function defined in the whole complex plane 
with branch points at z = 1 and z — oo. Instead, tech- 
niques which permit an analytic continuation from the 
interior to the exterior of the unit circle are needed. 

Sequence transformations can be used to accelerate the 
convergence of a hypergeometric series 2 Fi or to sum it in 
the case of divergence, which corresponds to an analytic 
continuation. Let us for example consider the following 
elementary special case of a hypergeometric function 2 Fi 
(Eq. (15.1.3) of M): 



ln(l + z) 



(-l) m z 



^-^ m + 1 

m=Q 

z 2 F 1 {l,l;2--z) 



(3.2) 



The infinite series converges only for \z\ < 1 , whereas the 
logarithm is with the exception of the cut along — 00 < 
z < — 1 defined in the whole complex plane. 
In Table 



Wynn's epsilon algo rithm, Eq. (2.5) 
Brezinski's theta algorithm, Eq. (A9), Levin's transfor- 
mation dfc (C, s n ), Eq. ( ]A22D , and the closely related se- 



quence transformation £>J^(C, s n ), Eq. ( |A25| ), are applied 
to the partial sums 



Sn(z) 



E 

771=0 



{-l) m z 



(3.3) 



of the hypergeometric series ( p.2[ ) for z = 7/2. The ap- 
proximations to the limit in Table | were always chosen 
in such a way that the transforms with the highest possi- 
ble transformation order were computed from a given set 
of input data. Thus, in the case of the epsilon algorithm, 
the appr oximations to the limit were chosen according 
to (2.12), in the case of th e theta algorithm, they were 
chosen according to ( A13 ), and in the case of the se- 
quence transformations rf j (C? s n) and 5i n ^(£, s n ), they 
were chosen according to ( Ajjj V. 

The second column of Table |, which displays the par- 
tial sums ( p. 3] ) , shows that the hypergeometric series 
(3.2) for ln(l + z) diverges quite strongly for z = 7/2. 
Nevertheless, it is apparently possible to sum this diver- 
gent series to its correct value. 

The results in Table | also show that Wynn's epsilon 
algorithm, which in the case of a p owe r series produces 
Pade approximants according to (2.7), is contrary to 
a widespread belief not necessarily the most powerful 



transformation. The theta algorithm and in particu- 
lar the two Levin-type transformations '(£, s n ) and 

^fe (Ci s ti)? which use the first term neglected in th e par- 
tial sum as a remainder estimate according to ( A21), pro- 
duce significantly better summation results. 

The other sequence transformations discussed in Ap- 
pendix |A] give better results than Wynn's epsilon algo- 
rithm, but are less effective than the transformations 

dk(Ci s n) and 5^ (C, s n ). For example, the two approxi- 
mations with the highest possible tra nsfo rmation orders, 
which Aitken's iterated A 2 process ( A4) produces from 
the partial sums sq(z), s±(z), sis(z) according to 
' ?|), are 



= 1.504 077 397 173. 
A$p = 1.504 077 396 169. 



(3-4) 
(3.5) 



Similarly, the iteration (A15) of B rezinski's theta algo- 
rithm produces according to ( Al9j ) the approximants 
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(2) 
i 

r(Q) 



1.504 077 404 830, 



J}"' = 1.504 077 394 094. 



(3.6) 
(3.7) 



These results show that sequence transformations can 
be very useful. Nevertheless, in the case of a real real ar- 
gument z it is actually not necessary to to use sequence 
transformations for doing the analytic continuations or 
for speeding up convergence. Instead, it is often simpler 
to exploit some known mathematical properties: Unless 
certain linear combinations of the parameters a, b, and c 
are positive or negative integers, a hypergeometric func- 
tion 2-F1 (a, b; c; z) can be expressed as the sum of two 
other 2-Fi's with a transformed argument w = 1 — z, 
w = l/z, w = 1/(1 — z), or w = 1 — 1/z, respectively. 
Thus, the argument w of the two resulting hypergeomet- 
ric series can normally be chosen in such a way that the 
two new series in w either converge, if the original se- 
ries in z diverges, or that they converge more rapidly if 
the original series converges too slowly to be numerically 
useful. 

For example, if |1 — z\ < 1 and if c — a — b is not a 
positive or negative integer, then we can use the analytic 
continuation formula (Eq. (15.3.6) of 



2 Fi(a,b; c;z) 
_ T(c)T(c -a-b) 
~ T(c-a)T(c-b) 
T(c)T(a + b- c) 



2 F 1 (a, b; a + b — c + 1; 1 — z) 

{l-z) c - a - b 



T(a)T(b) 

x 2 F 1 (c-a,c-b;c-a-b + l;l-z). (3.8) 

Let us now assume that z is only slightly smaller than 
1. Then, the convergence of the original hypergeometric 
series 2 Fi(a, b; c; z) will be very bad. However, the two 
hypergeometric series on the right-hand side with argu- 
ment 1 — z will converge rapidly in the vicinity of z = 1 . 
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With the help of this or similar analytic continuation 
formulas it is normally possible to compute a hyperge- 
ometric function 2 Fi(a,b; c; z) with real argument z ef- 
fectively since it is possible to find for every argument 
z € (—00, +00) two hypergeometric series with an argu- 
ment \w\ < 1/2 (see p. 127 of @ or Table I of ||). 

Unfortunately this approach does not necessarily work 
in the case of complex arguments z. Consider the points 



Zl,2 



1± V3 



(3.9) 



which both lie on the boundary of the circle of con- 
vergence because of |zi,2| — 1- in practice, it is either 
impossible or not feasible to evaluate a nonterminating 
hypergeometric series 2 F\ by adding up its terms if its 
argument z lies on the boundary of the unit circle. Con- 
sequently, something has to be done to speed up con- 
vergence or to accomplish a summation in the case of 
divergence. Unfortun atel y, the analytic continuation for- 
mulas of the type of (3.8) do not improve the situation if 
z = z\ 2 since 



1 - Z\,i 

1/(1-21,2) 
1 - lAi,2 



Z 2,l 
22,1 
21,2 
21,2 



(3.10a) 
(3.10b) 
(3.10c) 
(3.10d) 



Hence, the analytic continuation formulas cannot help if 
the z is close to zi, 2 - 

However, sequence transformations work for z = zi t 2- 
Let us consider the following elementary special case of 
a hypergeometric function 2 F\ (p. 38 of EJ): 



(1 + *)(!- 2) 



-2ce-l 



2 F 1 (2a,a + l;a;z). (3.11) 



For a > 0, the hypergeometric series converges in the 
interior of the unit circle, but it diverges on its boun dary. 
In Table ||, Wynn's epsilon algorithm, Eq. (2J3), and 



Levin's transformation ^[""'(C, s n ), Eq. ( A22), are used 



to sum the hypergeometric series (3.11) with a = 1/3 on 
the boundary of the unit circle, i.e., they are applied to 
the partial sums 



,(2) 



m=0 



(2/3) m (4/3) r 
(l/3) m m! 



(3.12) 



with z = zi = (1 + i\/3)/2. 

The partial sums ( [3.12] ) in the second column of Table 
[ij display a very unusual sign pattern and grow slowly in 
magnitude. Nevertheless, both the epsilon algorithm as 
well as the Levin transformati on ar e apparently able to 
sum the hypergeometric series (3.11) for z = (l + iv3)/2. 

The other sequence transformations discussed in Ap- 
pendix |A| are appa rently also able to sum the hypergeo- 
metric series ( 3.11 ) for z = (l + i\/3)/2. They give better 
results than Wynn's epsilon algorithm, but are less ef- 
fective than Levin's transformations djj, (Q, s„ ). The two 



approximations with the highest possible transformation 
orders, that can be produced from the partial sums s (z), 
s\(z), . . . , si5(z) by Aitken's iterated A 2 process (|A4|) , b y 
Brezinski's theta algorithm (A9) and it s iter ation ( |A15| ). 
and by the Levin-type transformation (A25), are 



(0) 



4 X) 



q(2) 



,(0) 
'10 



- 1.113 340 798 057 -I- i 1.326 827 896 288 . 

i) _ 

7 — 

- 1.113 340 798 408 + i 1.326 827 896 424 , 

= 

- 1.113 340 797 528 + i 1.326 827 893 689 . 
1.113 340 799 160 + i 1.326 827 895 539. 



3< 



(2) 



- 1.113 340 798 249 + i 1.326 827 894 967. 

3t = 

- 1.113 340 798 314 + i 1.326 827 895 955 . 

- 1.113 340 798 314 + i 1.326 827 895 955 . 

^(MoW) = 

- 1.113 340 798 414 + i 1.326 827 896 325 . 



(3.13) 
(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 



The numerical results presented here should suffice to 
support the claim of the author that sequence transfor- 
mations can be extremely useful numerical tools for the 
evaluation of special functions in general |l2|-[l7j and for 
the evaluation of Gaussian hypergeometric series 2 -Pi in 
special. Moreover, the numerical results shown above in- 
dicate that a computational algorithm, which would be 
capable of evaluating of a hypergeometric function 2^1 
with essentially arbitrary complex argument z and pa- 
rameters a, b, and c, should be a suitable comb inat ion of 
analytic continuation formulas of the type of (3^) with 
sequence transformations. 

Here it should be taken into account that a hyperge- 
ometric series 2^1 depends on four essentially arbitrary 
complex quantities a, b, c, and z. This makes the devel- 
opment of a general algorithm for its computation diffi- 
cult since a very large variety of special cases and com- 
putationally different situations have to be taken into 
account. Consequently, the development of such an al- 
gorithm would first require extensive numerical studies 
which - although undeniably interesting - would clearly 
be beyond the scope of this article. 



IV. 



HYPERGEOMETRIC SERIES WITH A 
NEGATIVE THIRD PARAMETER 



It is a direct consequence of its definition (3.1) that 



a hypergeometric series 2-^1(0, b; c; z) terminates after a 
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finite number of terms if either a or b is a negative integer. 
Moreover, 2-Fi(a, b; c; z) makes sense if both a and c are 
negative integers such that a — — m and c = —m—k with 
k, m = 1, 2, . . . : 



have c = —7/2. Thus, in T able IV we use Aitken's it- 
erate d A 2 process, Eq. (A4), Wynn's epsilon algorithm, 
Eq. (2.5), and Brezinski's theta algorithm, Eq. (A9), for 
the acceleration of the convergence of the partial sums 



2 F 1 (~m,b; 



-TO 



k;z) 



^ (-m- fc) M /x! 



(4.1) 



If c is a negative integer and if nei ther a nor b is a negative 
integer, then it follows from (3.1) that the hypergeomet- 
ric series 2F1 (a, b\ c; z) is in general undefined. However, 
the following limit exists for to = 0, 1, 2, . . . (p. 38 of 

lim ^vT2-Fi(a,&;c;z) = — —— 

c^-m 1 (cj [m + 1J! 

x 2 Fi(a + to+ 1,6 + to+ 1;to + 2; z) . (4.2) 



If c is not a negative integer, there are no problems 
with Pochhammer symbols in the denominators of the 
terms of the series that could become zero. Nevertheless, 
unpleasant numerical problems occur even if c is just a 
nonintegral negative real number. These problems can 
be demonstrated convincingly by trying to accelerate the 
convergence of two hypergeometric series 2F\{a, b; c; z) 
which differ only by the sign of the third parameter c. 

For that purpose, we consider in Table |l| and in Tables 
pV| and respectively, the hypergeometric series with 
a = 3/7, b = 5/2, z = 77/100, and c = ±7/2. 

In Table |l|, Wynn's epsilo n al gorithm, Eq. (2.5), 
Brezinski's theta algorithm, Eq. ( A9), and the Levin-type 



transformations % (Cj s n) 5 Eq. ( A22), and 6 y k rl '((,s ri 



An), 



Eq. ( |A25| ), are applied to the partial sums 



Sn(z) 



E 

m=0 



(3/7) m (5/2) 



(7/2), 



m m 



(4-3) 



of the hypergeometric series 2-Pi(3/7, 5/2; 7/2; z) with 
z = 77/100. The results in Table show that the con- 
vergence of this series can indeed be accelerated quite 
effectively by sequence transformations. This is also true 
for the other sequence transformations discussed in Ap- 
pendix |A|. The approximations with the highest trans- 
formation orders, that can be obtained from the partial 
sums sq(z), Si(z), siq(z) by Aitken 's iterated A 2 
process (A4) and by the iteration ( A15| ) of Brezinski's 
theta algorithm, are 



4 X) 



A 



(o) _ 

8 — 

(0) 



1.463 807 099 629 , 
1.463 807 099 563. 
1.463 807 143 254, 



jP = 1.463 807 103 421 . 



(4.4) 
(4.5) 
(4.6) 
(4.7) 



In Tables IV and |v|, we now consider a hypergeomet- 
ric series 2 F\ which is identical with the one in Table III 
except that its third parameter is negative, i.e., we now 



s„(z) 



n 

E 

m=0 



(3/7) m (5/2) t 
(-7/2) m m! 



(4.8) 



of the hypergeometric series 2-F1 (3/7, 5/2; —7/2; z) with 
z = 77/100, and in Table we use the iteration of Brezin- 
ski's theta algorithm, Eq. ( |A15 ), and the Levin-type 



tran sform ations d k (C,s n ), Eq. (|A22|). 



and 4 n) (C,s„), 

Eq. ( |A25|) . 

So far, the transformation results had always been very 
good as well as very reliable. In contrast, the results in 
Tables and are very bad and not reliable at all. 
For small transformation orders n, all transformations 
produce results which are by 5 orders of magnitude too 
small, and in the case of Brezinski's theta algorithm, the 
wrong results even seem to have converged with an ac- 
curacy of 4 decimal digits. For increasing transforma- 
tion orders n, there occur sudden and unmotivated sign 
changes, and only if n approaches 30, at least the epsilon 
algorithm in Table IV and the Levin transformation 
in Table ^ converge to the correct result. For n = 30, 
the other transformations show no indication of conver- 
gence and produce results which are still by some orders 
of magnitude too small. 

How can the disturbingly bad performance of sequence 
transformations in Tables IV and [v| be explained. The 
terms of a hypergeometric series 2-F1 (a, b; c; z) satisfy the 
2-term recursion 



n+l 



(a) n+1 (b) n+1 z 
(c) n+1 (n + l)\ 

_ (a + n)(b + n)z (a) n (b) n z r ' 



(c + n)(n + l) (c)„n! 



(4.9) 



Obviously, the factor (a+n)(b+n)z /[(c+n)(n+l)] on the 
right-hand side determines whether the terms increase 
or decrease in magnitude with increasing n. As long as 
this factor is greater than one in magnitude, the terms 
increase with increasing n, and as soon as this factor is 
smaller than one, the terms decrease. 

Thus, we only have to determine those values of n 
which satisfy 



(a + n)(b + n)z 



(c + n)(n + l) 



= 1 



(4.10) 



for given a, b, c, and z in order to find out for which 
values of the index n the terms of a hypergeometric series 
2-Fi(a, 6; c; z) change their growth pattern. 

In the case of the hypergeometric series in Table III 
with a = 3/7, b = 5/2, c = 7/2, and z = 77/100, there 
is no n > which satisfies condition ( 4.1 0| ) . This implies 
that the terms of this series decrease monotonously in 
magnitude with increasing n > 0. Moreover, the results 



7 



in Table II] show that the convergence of this series can 
be accelerated quite effectively. 

In the case of the hypergeometric series in Tables IV 
and with c = —7/2 , the situation is more complicated 
since condition (4.10) is satisfied by n w 22. Thus, the 
terms of this series initially increase up to n = 22, and 
only for n > 22 they decrease and ultimately pro duc e a 
convergent result. Accordingly, the partial sums (4.8) of 
this hypergeometric series initially look like the partial 
sums of a mildly divergent series, and only for n > 22, 
they behave like the partial sums of a convergent series. 
Therefore, it should not be too surprising that sequence 
transformations perfor m p oorly if they use as input data 
only the partial sums (4.8) with n < 22. However, even 
for 22 < n < 30, only the epsilon algorithm and the Levin 
transformation ultimately converge to the correct re- 
sult. This provides strong evidence that the irregular in- 
put data with small indices n have a detrimental effect on 
the performance of sequence transformations with large 
transformation orders, which also use input data with a 
correct behavior. 

We can test the hypothesis, that the irregular behavior 
of the initial partial sums (|l|) with n < 22 leads to the 
poor performance of sequence transformations in Tables 
[V and by skipping the terms up to n = 22 in the 
transformation processes. Accordingly, Wynn's epsilon 
algorithm, Eq. (|2.5|), and the Levin- type transformations 



VEq. ( ggp , and 4 (C,s„), Eq. ( ggp , use in 



in Table |Vj the modified partial sums 

(3/7) m (5/2) 



s (22; 



n+22 

E 



(-7/2) m m! 



m z m 



(4.11) 



of the hypergeometric series 2-Pi(3/7, 5/2; — 7/2; z) with 
2 = 77/100 as input data. This approach is possible since 
all sequence transformations c onsid ered in this article are 
quasi-linear, i.e., they satisfy (A30). 

The results in Table [v| indeed confirm our hypothesis 
since they are nearly as good as the results in Table [I], 
at least with respect to the transformation orders that 
are needed to achieve a given relative accuracy. 

The other sequence transformations discussed in Ap- 
pendix [X] produce results that are less good than those 
shown in Table VI. The two approximations with 



the highest possible transformation orders, that can be 



produced from the partial sums s^ 2 \z), sf^'iz) 



,( 22 )/ 



s^o'(z) by Aitken's iterated A 2 process ( A4 ) , by Brezin- 



ski's theta algorithm (A9) and its iteration (A15), are 



A (1) 


= 1.010 147 722 439 


10+ 5 , 


(4.12) 


-^10 


= 1.010 147 537 701 


10+ 5 . 


(4.13) 


W 12 


= 1.011 462 051 628 


10+ 5 , 


(4.14) 


fl(2) 
^12 


= 1.011 462 011 501 


10+ 5 , 


(4.15) 




= 1.010 233 908 825 


10+ 5 , 


(4.16) 


Jj 1} 


= 1.010 176 054 786 


10+ 5 . 


(4.17) 



The results in Table VI show that a hypergeometric 
series 2-P1 with a negative third parameter can be eval- 
uated reliably with the help of sequence transformations 
if the nonregular leading terms are excluded from the 
transformation processes. Unfortunately, this may lead 
to new problems since the number of terms, that initially 
grow in magnitude, may become quite large, in particu- 
lar if z is close to one. For example, if we increase the 
argument of the hypergeometric series in Table VI from 
z = 77/100 to z = 87/100 or to z = 97/100, then the 
number of terms, which initially grow in magnitude and 
have to be skipped, grow from n = 22 to n — 40 or even 
to n = 179. Moreover, more negative values of c also 
increase the number of terms that have to be skipped. 
For instance, if we consider the hypergeometric series 
2-F1 (3/7, 5/2; -13/2; z) with z = 77/100, z = 87/100, 
or z = 97/100, then we obtain n = 35, n — 63, or 
n = 279, respectively. These examples show that it is 
in principle possible to construct hypergeometric series 
2.F1 which can only be evaluated reliably with the help of 
sequence transformations if a very large number of terms 
is skipped in the transformation process. 

In the case of a Gaussian hypergeometric series, this 
poses no unsurmountable problems. Firstly, it is a triv- 
iality to compute the terms, even for very large indices. 
Consequently, a brute force evaluation of a Gaussian hy- 
pergeometric series by adding up the terms is possible as 
long as the argument is not too close to the boundary of 
the circle of convergence. Secondly, the highly developed 
mathematical theory of these functions makes it possible 
to simplify the numerical task with the help of known 
transformation formulas. For example, if the number of 
terms, that have to be skipped, is large because the ar- 
gument z of the hypergeometric series with a negative 
third parameter is close to one, then it may be a good 
idea to use the analytic continuation formula ( |3.8| ) . This 
would not necessarily solve the principal problems due 
to a negative third parameter, but the argument 1 — z of 
the two new hypergeometric series would then be small. 

However, the probably simplest approach would be the 
use of recurrence formulas. In this approach, one would 
have to evaluate two hypergeometric series 2 -Pi with suit- 
able positive values of the third parameter, for example 
with the help of sequence transformations, and to com- 
pute recursively the numerical value of desired hypergeo- 
metric series with a negative third parameter. Of course, 
many recurrence formulas are known. Nevertheless, nu- 
merous new three-term recurrence formulas satisfied by 
the Gaussian hypergeometric function 2 -Pi (a, b; c; z) are 
derived in Appendix |^. 

Unfortunately, these alternative approaches are in gen- 
eral not available if we have to evaluate an infinite series 
whose terms are determined numerically and only behave 
like the terms of a hypergeometric series 2 Pi (a, b; c; z) 
with a negative third parameter. In such a case, it can- 
not be excluded that the number of terms, that have to be 
skipped, would be so large that their computation would 
no longer be feasible. Then, neither the conventional 
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process of successively adding up the terms of the series 
nor sequence transformations would be able to provide 
reliable approximations to the value of such an infinite 
series. 



V. SUMMARY AND CONCLUSIONS 

A sequence transformation is a rule T which trans- 
forms a sequence {Snj^Lg with the (generalized) limit 
s to another sequence {s^j^o having the same (gen- 
eralized) limit s but different remainders r' n = s' n — s. 
The transformation process was successful if the new se- 
quence has better numerical properties than the origi- 
nal sequence. For example, a sequence transformation T 
accelerates convergence, if the transformed remainders 
r' n = s' n — s converge more rapidly t han the original re- 
mainders r n — s n — s according to (1.4), and T sums a 
divergent sequence, whose remainders r n do not vanish 
as n — > oo, to its generalized limit s if the transformed 



If the sets {s„, s„+i, . . . s„+z} of input data with increas- 
ing I have a sufficiently small starting index n, then all 

(n) 

transforms Ti will be affected by irregular input data, 



remainders r' n approach zero as n — > oo. 

As discussed in Section O, a sequence transformation 
T can be represented by an infinite set of doubly indexed 

(n) 

quantities TV; with k, n > that can be displayed in a 
two-dimensional array called the table of T . The super- 
script n denotes the minimal index occurring in the finite 
string {s n , s„ + i, . . . s n+ i} of sequence elements used for 
the computation of a given . The subscript k - usu- 
ally called the order of the transformation - is a measure 
for the complexity of the transformation process which 
yields T fc (n) . 

A convergence acceleration or summation process tries 
to obtain a better approximation to the (generalized) 
limit of the input sequence by proceeding on a certain 
path in the table of T. There is an in principle un- 
limited variety of different paths, but in practical ap- 
plications either order-constant path s defined in ( |2.3| ) or 
index-constant paths defined in (2.4) are normally used. 

On an order-constant path, a set {s n , s n +i, . . . s n +;} of 
input data of fixed length is used and the starting index 
ri of this set is increased successively. In contrast, an 
index-constant path uses sets of input data of increasing 
length, and it is always tried to compute from a given set 
of input data that element which has the highest 
possible transformation order k. 

Order-constant and index-constant paths differ sub- 
stantially. For example, on an index-constant path the 
available information is exploited more efficiently than 
on an order-constant path. Accordingly, index-constant 
paths normally produce better transformation results. 
Moreover, order-constant paths cannot be used for the 
summation of a divergent sequence since increasing n in 
the set {s„, s n +i, . . . s n +;} of input data normally only 
increases divergence. 

If, however, the leading elements of the input sequence 
{s„}^L behave irregularly, the principal advantages of 
index-constant paths can easily turn into disadvantages: 



albeit to a different degree. As shown by Tables IV and 
[v|, this can lead to unreliable or even completely nonsen- 
sical transformation results. In such a case, it is necessary 
to exclude the irregular input data from the transforma- 
tion process. This can be accomplished by using either 
an order-constant path or an index-constant path with a 
sufficiently large starting index. 

In Section III, the Gaussian hypergeometric function 
2-FT(a, b; c; z) is discussed, which is in general a multi- 
valued function defined in the whole complex plane with 
branch points at z = 1 and oo. However, it is defined by 



the power series (3.1) which only converges for \z\ < 1. 
Accordingly, sequence transformations can either be used 
for speeding up convergence or for accomplishing an ana- 
lytic continuation in the case of divergence. In Table |, an 
alternating hypergeometric series for ln(l + z) is summed 
effectively by sequence transformations for an argument 
z = 7/2 which is is far away from the unit circle, and 
in Table |l|, another special hypergeometric series, which 
converges only in the interior of the unit circle, is evalu- 
ated for an argument z = (1 + iV3)/2 that is located on 
the boundary of the unit chicle. 

For most parameters a, b, and c, sequence transforma- 
tions greatly facilitate the evaluation of a hypergeometric 
series 2-F1, and the good transformation results presented 
in Tables | and [n] are fairly typical. However, as discussed 
in Section there is an important and instructive ex- 
ception: If the third parameter c of a hypergeometric 
series is a negative real number, then the terms of this 
series initially increase in magnitude and look even for 
\z\ < 1 like the terms of a mildly divergent series. Only 
for sufficiently large values of the index, the terms de- 
crease and ultimately produce a convergent result. 

The use of these irregular terms as input data seri- 
ously affects the performance of sequence transforma- 
tions and leads to unreliable and sometimes even com- 
pletely nonsensical results. This is demonstrated by ap- 
plying in Table III and in Tables IV and respectively, 
sequence transformations to the partial sums of the hy- 
pergeometric series 2-^1(0, b; c; z) with a = 3/7, b = 5/2, 
z = 77/100, and c = ±7/2. 

In Table III, sequence transformations are applied to 



the hypergeometric series with the positive value of the 
third parameter. As expected, the transformation results 
are very good. However, in Tables IV and [v|, where the 
hypergeometric series with the negative value of the third 
parameter is considered, the transformation results are 
both unreliable and bad. 

Nevertheless, it is possible to compute the hypergeo- 
metric series 2 F 1 (3/7, 5/2; -7/2; 77/100) efficiently and 
reliably with the help of sequence transformations. How- 
ever, one cannot use an index-constant path with a small 
minimal index, as it was done in Tables IV and [v| In- 
stead, one should either use an order-constant path or an 
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index-constant path with a sufficiently large minimal in- 
dex, as it was done in Table VI , where all irregular terms 
were excluded from the transformation processes. 

The sequence transformations, which are used in this 
article, are all described in Appendix |A|. The use of 
several transformations was quite intentional. The au- 
thor wanted to make clear that problems due to irregu- 
lar input data are not restricted to some special sequence 
transformations only. Of course, the results in Tables IV 
and ^ show that different sequence transformations re- 
spond differently to irregular input data. However, this 
is quite helpful and can protect us against misinterpre- 
tations. For example, in Table IV Brezinski's theta al- 



gorithm produced transformation results which seemed 
to have converged with an accuracy of 4 decimal digits, 
but were actually by 5 orders of magnitude too small. 
Fortunately, the other transformations in Tables [V and 
^ produced different results, which provided strong ev- 
idence that the transformation results were unreliable. 
Consequently, it is recommendable to use in convergence 
acceleration and summation processes more than a sin- 
gle transformation whenever possible. This is particu- 
larly important if numerically determined data are to be 
transformed, about which very little is known. If sev- 
eral different sequence transformations produce consis- 
tent results, then it is very likely that these results are 
indeed correct although it is of course clear that purely 
numerical results cannot be a substitute for a rigorous 
mathematical proof. 

It looks like a contradiction t hat in Appendix [§] only 
the Levin-type transformation (A25) was used for the 
computation of the infinite coupling limit of the sex- 
tic anharmoni c os cillator. However, the divergent pertur- 
bation series (]B9|), whose leading coefficients also show 
an irregular behavior, constitutes a very demanding sum- 
mation problem, for which the other transformations dis- 
cussed in Appendix ^ are not powerful enough. 

Here, it should be emphasized once more that the two 
examples considered in this article - the Gaussian hy- 
pergeometric series 2-Fi (a, b;c;z) with a negati ve t hird 
parameter and the divergent perturbation series ( |B9| ) for 
the infinite coupling limit - are comparatively sim- 
ple model problems, and the leading irregular terms of 
their series expansions pose no unsurmountable compu- 
tational problems. This is largely due to the fact that 
it is relatively easy to find out which terms beha ve ir - 
regularly. In the case of the perturbation series (p39|), 
we are in the fortunate situation that the leading large-n 

,(3) 



asymptotics (BIO) of the coefficients c„ is known, and 
in the ca se of the hypergeometric series, one only has 
to solve ( 4.1 0| ) in order to find out for which indices n 
the terms change their growth pattern. Moreover, in the 
case of a 2-F1 there are numerous alternative computa- 
tional approaches. For example, with the help of recur- 
rence formulas the evaluation of Gaussian hypergeomet- 
ric series with a negative third parameter can be avoided 
completely. Consequently, in Appendix ^ numerous new 



recurrence formulas satisfied by a 2-F1 (a, 6; c; z) are de- 
rived. 

Finally, the author wishes to express his hope that this 
article will inspire additional research on the evaluation 
of special functions with the help of sequence transfor- 
mations. There can be no doubt that sequence trans- 
formations are normally extremely useful tools for the 
evaluation of special functions, and since the terms of 
the series expansions for special functions are compara- 
tively simple and explicitly known, we can even hope to 
gain additional insight from those cases in which sequence 
transformations fail to produce good transformation re- 
sults. 
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APPENDIX A: SEQUENCE 
TRANSFORMATIONS 

This appendix gives a short description of all the se- 
quence transformations that are used in this article. Fur- 
ther details plus additional references can be found in 
P,RJT3]. Here, the same notation as in |E| is used. 

One of the oldest sequence transformations (see for in- 
stance pp. 90 - 91 of [§) is Aitken's A 2 formula (48): 



4 n) 



[As T 



A 2 i 



(Al) 



The (forward) difference operator A acts for all integers 
n > on a function f{n) according to 



A/(n) = /( n + !)-/(«). 



(A2) 



The Aitken formula (|Al[ ) is by construction exact for 
model sequences of the type s n = s + cA™ with c ^ 
and A ^ 1. If the numerical values of three consecutive 
elements s„, s n +i, and s n+ 2 of this model sequence are 
known, then the (generalized) limit s of this sequence can 
be computed according to 



(n) 



(A3) 



no matter whether the sequence converges (|A| < 1) or 
diverges (|A| > 1). 

The power and practical usefulness of Aitken's A 2 for- 
mula is of course limited since it is designed to eliminate 
only a single exponential term from the elements of the 
model sequence mentioned above. However, the quanti- 



ties can again be used as input data in ( Al). Hence, 
the A 2 process can be iterated, yielding the following 
nonlinear recursive scheme fll, Eq. (5.1-15)]: 
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(n) 



A 



A (n) - A (n) 



A 2 ^' 



(n) 



(A4a) 
(A4b) 



In this article, the difference operator A acts only on the 
superscript n and not on the subscript k of a doubly in- 
dexed quantity like Au n . A more detailed discussion of 
Aitken's iterated A 2 process as well as additional refer- 
ences can for instance be found in Section 5 of [|l2| or in 
p9| . The iteration of other sequence transformations is 
discussed in Q. 

In the case of Aitken's iterated A 2 process, the ap- 
proximation to the limit of the input sequence with the 
highest possible transformation order depends upon the 
index to of the last sequence element s m which was used 
in the recursion. If to is either even or odd, m = 2/i or 
m = 2/i + 1, respectively, the approximations to the limit 
are chosen according to 



{«o, si, 

{Sl, 82, ■ ■ 



> s 2fi+l} 



A {0) 

A 1 

A w 



(A5) 
(A6) 



As in the case of Wynn's e algorithm, these two rela- 
tionships can with the help of the notation \x\ for the 
integral part of x be expressed by a single equation (Eq. 
(5.2-6) of @): 



2[[m/2]]) 5 m _2[m/2]] + lj ■ 



A 



(m-2[m/2]) 
ffm/2] 



(A7) 



The behavior of many practically relevant convergent 
sequences {s n }%L can be characterized by the asymp- 
totic condition 



lim 



Sn+l - S 



f). 



(A8) 



where s = is the limit of the sequence {snj^Q. This 
condition closely resembles the well known ratio test for 
infinite series. A convergent sequence satisfying (AS) 
with \p\ < 1 is called linearly convergent, and it is called 
logarithmically convergent if p = 1 . 

It is one of the major weaknesses of the other wise 
very powerful and very useful epsilon algorithm ( |2.5| ) 
that it does not work in the case of logarithmic conver- 
gence. Brezinski showed that this principal weakness can 
be overcome by a suitable modification of the recursive 
scheme ( |2.5| ), which leads to the so-called theta algorithm 



(n) 

o(") 
7 2fc+l 

n(n) 
7 2fc+2 



= o, 



q(n) 
7 



g(«+l) 
7 2fc-l 

fl("+l) 
7 2fe 



[a^: +i) ] 



A 2 e». 



(A9a) 
(A9b) 

(A9c) 



As in the case of Aitken's iterated A 2 process (A4), it is 
assumed that the difference operator A acts only upon 
the superscript n and not on the subscript k. 

Again, the approximation to the limit of the input se- 
quence depends upon the index to of the last sequence 
element s m which was used in the recursion. If we have 
to = 3/i, to = 3/1+1, or to = 3/i + 2, respectively, the 
approximations to the limit with the highest transforma- 
tion orders are chosen according to 



{so,si, 
{si,s 2 , ■ • 

{s 2 , S3, ■ ■ 



•-,S 3m } 
. s 3ai+i} 
> S3 M +2J 



q(0) 
7 2 M 

7 2/i 

g(2) 
7 2n 



(A10) 
(All) 
(A12) 



These three relationships can be expressed by a single 
equation (Eq. (10.2-8) of ||): 



{ s m-3[m/3]i s m-3[m/3] + l; 



q(m-3[m/31) 
2[m/3] 



(A13) 



Further details on the theta algorithm as well as addi- 
tional references can be found in Section 2.9 of Q or in 
Sections 10 and 11 of @. 

As for example discussed in |5(i|| , new sequence trans- 
formations can be constructed by iterating explicit ex- 
pressions for sequence transformations with low trans- 
formation orders. The best known example of such an 
iterated sequence tr ansf ormation is probably Aitken's it- 
erated A 2 process (AA) which is obtained by iterating 



(Al) 



The same approach is also possible in the case of the 
theta algorithm (|A9|). A suitable closed- form expression, 
which may be iterated, is (Eq. (10.3-1) of O) 



t) 



(n) 



Sn+l 

[As n ] [As n+1 ] [A 2 s n+1 ] 
[As„ +2 ] [A 2 s„] - [As n ] [A 2 s n+1 ] 



(A14) 



The iteration of this expression yields the following non- 
linear recursive scheme (Eq. (10.3-6) of |0|): 



<t(«) 

Jo 

J k+l 



= Jl 



(n+1) 



[aj^] [Ajr i} ] [a 2 jy +1> ] 



0+i)i 



r(«+l)l 



[Ajt +2) ] [A*jn - [Ajn [A^r 11 ] 



(Al5a) 



(A15b) 



2k+l 



In convergence acceleration and summation processes, 

(n) 

the iterated transformation J7" fe has similar properties 
as Brczinski's theta algorithm from which it was derived: 
Both transformations are very powerful as well as very 
versatile. J^ 1 is not only an effective accelerator for lin- 
ear convergence as well as able to sum divergent alternat- 
ing series, but it is also able to accelerate the convergence 
of many logarithmically convergent sequences and series 
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Again, the approximation to the limit of the input se- 
quence depends upon the index to of the last sequence 
element s m which was used in the recursion. If we have 
to = 3/i, m = 3/i + 1, or to = 3/i + 2, respectively, the 
approximations to the limit with the highest transforma- 
tion orders are chosen according to 



{s ,si, . . . ,s 3jU } 
{si,s 2 , • • . ,s 3jU+ i} 
{s 2 ,s 3 , . . . ,s 3ai+2 } 



-7(0) 
-7(1) 

Ji 2) ■ 



(A16) 
(A17) 
(A18) 



These three relationships can be expressed by a single 
equation (Eq. (10.4-7) of pf): 



m-3lm/3} 



s m-3[m/3]] + l! • 
(m-3[m/3l) 



[m/3] 



(A19) 



So far, only sequence transformations were considered 
which use as input data the elements of the sequence to 
be transformed. However, in some cases structural infor- 
mation on the dependence of the remainders r n on the 
index n is available. For example, it is well known that 
the truncation error of a convergent series with strictly al- 
ternating and monotonously decreasing terms is bounded 
in magnitude by the first term not included in the par- 
tial sum and has the same sign as this term (see for in- 
stance p. 132 of |2j). The first term neglected is also 
the best simple estimate for the truncation error of a 
strictly alternating nontcrminating hypergcometric series 
2 F (a, 0; -x) with a,P,x > (Theorem 5.12-5 of @). 
Such an information on the behavior of the truncation 
errors can be extremely helpful in a convergence acceler- 
ation or summation process. Unfortunately, the sequence 
transformations considered so far are not able to benefit 
from it. 

A convenient way of incorporating such an informa- 
tion into the transformation process consists in the use 
of remainder estimates {Lu n }^ =Q . Because of the explicit 
incorporation of the information contained in the remain- 
der estimates, sequence transformations of that kind are 
potentially very powerful and as well as very versatile. 

The best-known example of such a sequence transfor- 
mation is Levin's transformation |p8[ w hich is both very 
versatile and very powerful ||[l2 59^6lf|: 



3=0 



fc-1 



jj {C + n + k) 



k-l 



EM) 



3=0 



(C + n + jf- 1 
jj (C + n + k) k ~ 1 



'->n+j 



(A20) 



The shift para meter £ has to be positive in order to admit 
n = in ( A20). The most obvious choice, which is always 
used in this article, is C — 1- Recurrence formulas for the 



numerator and denominator sums of jC^ (£, s n , w n ) can 
be found in Section 7.2 of @. 

Levin's transformation is based on the assumption that 
the remainders r n of the input sequence can for all n > 
be approximated by a remainder estimate u> n , which 
should be chosen such that s n — s = ui n \c + 0(n -1 )] 
as tl — ► oo , multiplied by a polynomial in l/(n + £) with 
£ > 0. Levin p8[ introduced several simple remainder 
estimates for infinite series which give rise to several vari- 
ants of Levin's sequence transformation. Further details 
on Levin's transformation can for instance be found in 
Section 7 of @. 

In this article, we only consider the remainder estimate 



As r , 



ln+1 , 



(A21) 



which was first proposed by Smith and Ford j5Jj] • It yields 
the following variant of Levin's transformation (Eq. (7.3- 

9) of my. 



4 n) (C,s„) = £nC,s„,A s „) 



(«)/ 



(A22) 



Levin's transformation is based on the implicit as- 
sumption that the ratio [s n — s]/uj n can be expressed as 
a power series in 1/ (n + C) . A different class of sequence 
transformations can be derived by assuming that the ra- 
tio [s„ — s]/u) n can be expressed as a so-called factorial 
series according to (Section 8 of p2^|) 



= s + u n ^Cj/in + Oj 



(A23) 



3=0 



where (n + C)j = T(n + j + Q /T(n + £) is a Pochhammer 
symbol (p. 3 of §U). In this way, the sequence transfor- 
mation (Eq. (8.2-7) of @) 



(_1V ( k \ (C + » + j)fc-i s n+j 



3=0 



E 



3=0 



k\ (C + n + j) fc -i 1 
jj (C + n + k)k-i ojn+j 



(A24) 



can be derived which is formally very similar to Levin's 
sequence transformation. The only difference between 
the transformations C^\(,, s n ,uj n ) an d (C, s ni L0 n ) is 
that the powers (C + n + j) k ~ l in (A2C) are re place d 
by the Pochhammer symbols (C + n + j)k-i in ( A24| ) . 
Again, the shift parameter £ has to be positive in order 
to admit n = in (A24), and the most obvious choice 
is also C = 1 which is exclusively used in this article. 
Recurrence formulas for the numerator and denominator 
sums of (C, s. 



can be found i n Sec tion 8.3 of flL2\ 
If we use the remainder estimate ( A21 ) also in (|A24|) 



we obtain the following sequence transformation (Eq. 
(8.4-4) of 
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4 n) (C,Sn,A S „). 



(A25) 



It was shown in several articles that the transformation 



(|A24 ) as well as its variant ( A25 ) can be very effective 
||^,|l]j23f^,||-|l|,|2|-||, in particular if strongly di- 
vergent alternating series are to be summed. 

In the case of the transformations ( A22[ ) and ( A25| ), 
the approximation to the limit with the highest transfor- 
mation order is given by 



{so, si, 

7» 



l} 



(A26) 



where E\! l> ((, s n ) stands for either dfc™^(C,Sn) or 

4 n) (c, 5 „). 

If the input data s n are the partial sums of a (for- 
mal) power series for some function f(z) accordi ng to 
( |2.6|) , s n — fn(z), then the transformations ( A22 ) 
and ( |A25 ) produce rational functions d^f 1 (£, f n (z)) and 

^i™*' {(> fn(zj) i whose numerator and denominator poly- 
nomials are of degrees k + n and k in z, respectively (Eqs. 
(4.25) and (4.26) of pa]). Moreover, the rational approx- 



imants (C, f n (z)) and 5j™ ((,f n (z)) satisfy the fol- 
lowing asymptotic error estimates as z — > (Eqs. (4.28) 
and (4.29) of @), 



f(z)-d^(C,f n (z)) = 0(z fc +"+ 2 ), 
f(z) - *£°(C, /»(*)) = 0(z fe +"+ 2 ), 



(A27) 
(A28) 



which are very similar to the well known accuracy- 
through-order relationships satisfied by Pade approxi- 
mants |34| ]. 

It is a typical feature of all sequence transformations 
discussed in this Appendix that they are both homoge- 
neous and translative: If the elements of two sequences 
{s n }n=o and Wnj^Lo satisfy 

cr„ = as n + b , (A29) 

where a and b are suitable constants, then 

T(a n ,a n+ x,...) = aT(s n ,s n+1 , . . .) + b. (A30) 

Sequence transformations T satisfying this condition are 
called quasi-linear in the book by Brezinski and Redivo 
Zaglia ||. In Section 1.4 of this book, a detailed discus- 
sion of the properties of quasi-linear sequence transfor- 
mations as well as further references can be found. 



APPENDIX B: THE INFINITE COUPLING 
LIMIT OF THE SEXTIC ANHARMONIC 
OSCILLATOR 

The detrimental effect of the irregular behavior of the 
leading elements of a sequence in convergence accelera- 
tion and summation processes is not restricted to math- 
ematical model problems but occurs also in the mathe- 
matical treatment of scientific problems. 



This will be shown by performing extensive summation 
calculations for the so-called strong coupling limit £3 of 
the sextic anharmonic oscillator. The quartic (m = 2), 
sextic (m = 3), and octic (m = 4) anharmonic oscillators 
are defined by the Hamiltonians 



H(P) 



P 



fix 



2 m 



m 



2,3,4, (Bl) 



and the strong coupling limit k m of the ground state 
energy eigenvalue E^ m {(3) of this Hamiltonian is defined 
by 



lim fiW^/M) 



(B2) 



Ever since the seminal work of Bender and Wu [|66| 
the divergent weak coupling perturbation expansion 



£ (m) (/3) 



(B3) 



for the ground state energy of an anharmonic oscilla- 
tor has been considered to be the model example of a 
strongly divergent quantum mechanical perturbation ex- 
pansion which has to be summed in order to produce 
numerically useful results. Accordingly, there is an ex- 
tensive literature on the summation of the divergent per- 
turbation expansions of the anharmonic oscillators (see 
for example [T^] and references therein). 

In addition to the divergent weak coupling expansion 



(B3), there is also a strong coupling expansion 



n=Q 



It can be shown that this expansion converges for suffi- 
ciently large values of f3 |69|-|7l]]. However, the computa- 
tion of the perturbative coefficients Kn is very difficult 
(see for example p9|1 a n d ref erences th erein) . 

It follows from JB2|), (B3), and (B4) that the infinite 
coupling limit k m corresponds to the leading coefficient 
of the strong coupling expansion (B4) according to 



K, 



(,„) 



(B5) 



The weak coupling perturbation expansion ( |B3| ) can- 
not be used in a straightforward way for a calculation 
of the strong coupling limit k m . However, this can be 
accomplished comparatively easily with the help of the 
following renormalized weak coupling perturbation ex- 
pansion (Eqs. (3.30) - (3.31) of §§): 



E {m \(3) = (1-k) 



- 1/2 E 

n=0 



cl m) K B 



(B6) 



This expansion is based on a renormalization scheme in- 
troduced by Vinette and Cizek j72|. In this approach, 
the original coupling constant j3 €[0, 00) is transformed 
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into a renormalized and explicitly m-dependent coupling 
constant k S [0, 1) according to (Eq. (3.19) of (2fj) 

where (Eq. (3.17) of @) 

m(2m- 1)!! 



5,, 



(B8) 



For the sextic (to = 3) case, these expressions correspond 
to B 3 = 45/4 and [3 = 4k/ [45(1 - k) 2 ]. 

Thus, the infinite coupling limits k 3 of the sextic an- 
harmonic oscillator can be exp ressed by the renormalized 
weak coupling expansion ( pif ) according to (Eqs. (3.43) 
and (3.44) of ||) 



fc 3 = [45/4] 



OO 

1/4 E - 

n=0 



(3) 



(B9) 



The summation of either this or the perturbation series 
(B6) with m = 3, from which (B9) was derived, is a 
formidable computational problem. This follows at once 
from the large-n asymptotics of the renormalized per- 
turbative coefficients for the sextic anharmonic oscillator 
(Eq. (3.34) of ||): 



(128) 



1/2 



(-1) T 

VI ~ 

x T(2n+ 1/2) (64/[45tt 2 ])' 



(BIO) 



It should be note d that the summation of the pertur- 
bation series (B9) for k 3 is much more demanding than 
the summation of the divergent asymptotic expansions 
for special functions since their coefficients c n grow es- 
sentially like nl @-||,|6|. Although Pade approximants 
- or equivalently Wynn's epsilon algorithm - are in prin- 
ciple capable of summing alternating divergent power se- 
ries whose coefficients c„ grow essentially like (2n)! in 
magnitude, the convergence of Pade approximants is too 
slow to be practically useful. Moreo ver, i t was shown 
in [^H that the Levin transformation ( A22| ) produces in 
the case of the perturbation expansions for the anhar- 
monic oscillators sequences of approximants which ini- 
tially seem to converge but which ultimat ely d iverge. In 
contrast, the Levin-type transformation (A25) produces 
comparatively good results. 

Th us, in analogy to 25 we sum the perturbation se- 
ries ( p9| ) with the help of the Levin-type transforma- 
tion SjT (£, s n ) defined in (A25). It was shown in sev- 



eral articles that the sequenc e tra nsformation (A25) as 
well as the transformation (A24), from which it was 
derived, is apparently very effective, in particular if 
stro ngly d ivergent alternating series are to be summed 

In our summation calculations for k 3 we use the renor- 

(3) 

malized coefficients cj> with < v < 300. The coef- 
ficients were calculated using the exact rational arith- 
metics of Maple by solving a system of nonlinear dif- 
ference equations as described in the Appendix of p5[. 



Unfortunately, Eq. (A22) in the Appendix of |25j] , which 
specifies the system of nonlinear equations, contains a 
typographical error. Correct is 



4jG 



(n) 



2(j + 1)(2j + 1)G;.; ) 1 + 



1 



-G 



(n-l) 



G 



(n-l) 
j-1 



2 ]T G< fc) G 



(fe) n {n-k) 



(Bll) 



fe=i 



The topic of this article is the study of the impact of ir- 
regular input data on the performance of sequence trans- 
formations. Accordingly, we have to investigate whether 
the renormalized coefficients c n behave irregularly for 
small indices n. For that purpose, we list in Table VII 

(3) 

selected renormalized coefficients c„ as well as the cor- 
responding ratios 



Ci 3 > 



(-1) 



n+1 _2 (3) 
7T Cn 



457i- 



/ 128r(2n+ 1/2) V 64 



(B12) 



that are obtained by dividing the renormalized coeffi- 

(3) 

cients c n by the leading order of their large-n asymp- 



totics according to (B1C) 



The last column in Tabic VII shows quite clearly that 



(3) 

the renormalized coefficients c„ deviate for small val- 
ues of n considerably from their larger-order behavior. 

Firstly, the coefficients Cq 3 ^ and q apparently possess 

"wrong" sign. Secondly, the coefficients c n initially de- 
crease in magnitude, and only for n > 4 they grow as 
they should according to ( B1C ). Nevertheless, it is only 
a relatively mild irregularity, which affects only a few of 

(3) 

the available coefficients c v with < v < 300. 

(3) 

The impact of the irregular coefficients c n with small 
indices n can be checked by computing for I = 0, 1, 2, . . . 
and for n < 299 — I the approximants 



kt' 1 " = ^(1,8^) 



(B13) 



to the infinite coupling limit k 3 of the sextic anharmonic 
oscill ator. Here, Sn^ is the Levin-type transformation 



(A25), and 



,('1 



[45/4] 



n+l 

1/4 E 

u=0 



.(3) 



(B14) 



is a partial sum of the perturbation series (BE) which 
skips the first I terms in the transformation process. 



In Table [Villi, the approximants kJ with the three 



highest possible values of n < 299 — Z are listed for I < 12. 

If we compare the results with the extremely accurate 
result of Vinette and Cizek (Eq. (69) of @) ' 

k 3 = 1.144 802 453 797 052 763 765 457 534 149 549 , 

(B15) 
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which was obtained nonperturbatively, we see that we 
gain 5 decimal digit s by skipping the first 7 terms of the 
perturbation series (B9) for k%. For I > 8, the accuracy 
of the summation results deteriorates again. 

This is a very remarkable gain of accuracy, if we take 
into acco unt that the summation of the perturbation se- 
ries ( p9| ) is a formidable task and that the leading terms 

(3) 

Cn disp lay a n only relatively mild irregularity, as shown 
in Table VII. Nevertheless, the results show quite clearly 
that the transformation order is not the only criterion 
which affects the performance of a sequence transforma- 
tion. The results in Table VIII show that may be more 



effective to use smaller sets of more regular input data. 

APPENDIX C: RECURRENCE FORMULAS FOR 
THE GAUSSIAN HYPERGEOMETRIC 
FUNCTION 

Many three-term recurrence formulas satisfied by 
the Gaussian hypergeometric function 2*1 (a, 6; c; z) are 
known. For example, on pp. 557 - 558 of |39fl or on pp. 
46 - 47 of pi] , the following formulas can be found: 

(c - a) 2 *i(a - 1, 6; c; z) 

+ [2a — c — (a — b)z] 2*1 (a, b; c; z) 

+ a(z-l) 2 F 1 (a + l,b;c;z) = 0, (CI) 
(c - b) 2*1 (a, b - 1; c; z) 
+ [2b — c — (b — a)z] 2*1 (a, b; c; z) 

+ b{z-l) 2 F 1 (a,b+l;c;z) = 0, (C2) 
c(l - c)(l - z) 2*1 (a, b; c - 1; z) 
+ c[c — 1 — (2c — a — b — l)z] 2F1 (a, b; c; z) 

+ {c-a){c-b)z 2 Fi{a,b;c+l;z) = 0, (C3) 
(b - a) 2*1 (a, 6; c; z) + a 2 F 1 (a + l,b;c;z) 

- b 2 F 1 (a,b+l;c;z) = 0, (C4) 
(6 -c) 2 *i (a, 6- l;c; z) 
+ (c - a - b) 2 Fi (a, 6; c; z) 

+ a(l-z) 2 *i(a+l,6;c;z) = 0, (C5) 
c[a — (c — b)z] 2 Fi (a, 6; c; z) 
— ac(l — z) 2-F\(a + 1, c; z) 

+ (c- a)(c- b)z 2 F 1 (a,b-c+ l;z) = 0, (C6) 
(l-c) 2 F 1 (a,&;c- l;z) 
+ (c- a - 1) 2^1(0, b;c;z) 

+ a 2 F 1 (a+ l,6;c;z) = 0, (C7) 
(a - c) 2-Fi(a - 1, b; c; z) 
+ (c - a - 6) 2-F1 (a, 6; c; z) 

+ 6(1 - z) 2*1(0, 6 + 1; c; z) = , (C8) 
(a - c) 2 Fi (a - 1, 6; c; z) 
+ (c - b) 2 Fi(a, b - 1; c; z) 

+ (6-a)(l-z) 2 Fi(a,6;c;z) = 0, (C9) 
(-c) 2 Fi(o- 1,6; c;z) 



+ c(l - z) 2*1 (a, b; c; z) 

+ (c-6)z 2-F1 (a, 6; c+l;z) = 0, (C10) 
(c - a) 2*1 (a - 1, 6; c; z) 

- (c-l)(l-z) 2 *i(a,6;c-l;z) 

+ [a - 1 - (c - 6 - l)z] 2 Fi (a, 6; c; z) = , (Cll) 
c[b - (c - a)z] 2*1 (a, 6; c; z) 

- 6c(l — z) 2*1 (a, b + 1; c; z) 

+ (c- a)(c- 6)z 2*1 (a, 6; c+ l;z) = 0, (C12) 
(1 - c) 2 F 1 (a,b;c - l;z) 
+ (c - b - 1) 2*1 (a, 6; c; z) 

+ 62*1 (a, 6+ l;c;z) = 0, (C13) 
(-c) 2*1 (a, b- l;c;z) 
+ c(l - z) 2*1 (a, 6; c; z) 

+ (c- a)z 2*1 (a, 6; c+l;z) = 0, (C14) 
(c- 6) 2*1 (a, 6- l;c;z) 

- (c-l)(l-z) 2 *i(a,6;c-l;z0 

+ [b - 1 - (c - a - l)z] 2*1 (a, 6; c; z) = . (C15) 

It is a typical feature of these recurrence formulas that 
there is a hypergeometric function 2*1(0, 6; c; z) plus two 
other 2*1 's which differ with respect to only one of the 
three parameters by ±1. However, recurrence formulas, 
which contain a hypergeometric function 2 Fi(a, b; c; z) 
plus two other 2*1 's which differ with respect to two or 
even three parameters by ±1, can be constructed com- 
paratively easily. For that purpose, we combine those of 
the recurrence formulas given above, whose the third pa- 
rameter c assumes at least two different values, with the 
linear transformation formulas (see for example p. 559 of 
f§ or p. 47 of 0) 

2*1 (a, b- c; z) = (1 - zf- a - b 2*1 (c - a, c - b; c; z) (C16) 
= (1 - z)- a 2*1 (a, c - b; c; z/(z - 1)) (C17) 
= (1 - z)- b 2*1 (c - a, b; c; z/(z - 1)) . (C18) 

For th e de rivation of new recurrence formulas, we re- 
place in flC3p a by c — a and 6 by c — 6. This yields: 

c(l — c)(l — z) 2*1 (c — a, c — 6; c — 1; z) 
+ c[c — 1 — (a + 6 — l)z] 2*1 (c — a, c — 6; c; z) 
+ a6z 2*1 (c- a, c- 6; c+ l;z) = 0. (C19) 

If we now combine this relationship with the linear trans- 
formation ( C16| ), we obtain the following recurrence for- 
mula, where all three parameters of the hypergeometric 
functions change simultaneously: 

c(l-c) 2 *i(a-l,6-l;c-l;z) 
+ c[c — 1 — (a + b - l)z] 2*1 (a, b; c; z) 

+ abz(l- z) 2 Fi(a + l,b+l;c+l;z) = 0. (C20) 



If we now proceed in @, @), and in jC10| ) - ( |C15| ) in 
exactly the same way, we obtain the following recurrence 
formulas: 
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c(c — a — bz) 2-P1 (a, 6; c; z) 

- c(c - a) 2-P\(a - l,6;c; z) 

+ a6z(l-z) 2 Pi(a + l,6+l;c+l;z) = 0, (C21) 
(1 -c) 2 Fi(o- 1,6- l;c- l;z) 
+ (a - 1)(1 - z) 2 Fi(a, b;c;z) 

+ (c-a) 2 F 1 (a-l,b;c;z) = 0, (C22) 
c 2 Fi(a + 1, 6; c; z) 

- c 2 Pi(a, 6; c; z) 

-bz 2 F 1 (a + l,b + l;c+l;z) = 0, (C23) 

a(l — z) 2-Fi( a + 1) &j c 'i z ) 

+ (l-c) 2 Pi(a-l,6-l;c-l;z) 

+ [c-a-1- (b-l)z] 2 F 1 (a,b;c;z) = 0, (C24) 
c(c — 6 — az) 2-F1 (a, 6; c; z) 

- c(c - b) 2 Fi(a, b - 1; c; z) 

+ a6z(l - z) 2 Pi(a + 1,6+ l;c+ l;z) = 0, (C25) 
(l-c) 2 Fi(a-l,b- l;c- l;z) 
+ {b- l)(l-z) 2 Fi(a,6;c;z) 

+ (c-6) 2 Fi(o,6- l;c;z) = 0, (C26) 
c 2 Fi(a, b + 1; c; z) 

- c 2 Fi(a,b;c; z) 

- az 2 F 1 {a+ 1,6+ l;c+ l;z) = 0, (C27) 
6(1 - z) 2 Fi(a,6+ l;c;z) 
+ (l-c) 2 Pi(a-l,6-l;c-l;z) 

+ [c - 6 - 1 - (a - l)z] 2 -Fi(a, 6; c; z) = . (C28) 

Not all of these recur rence formulas are independent. For 
example, ( C21 ) and ( C25 ) can be transformed into each 
other by i nterc hangi ng a a nd 6 . This is a lso tr ue fo r 
the pairs (|C2|) and (|C2|), (|c|f) and (|C27|), and dC2j) 
and (C2S), which can be transformed into each other by 



interchanging a and 6. 

For the derivation of recurrence formulas which d iffe r 
with respect to two parameters by ±1, we replace in ( |C3| ) 
a by c — a and z by z/(z — 1). This yields: 



c(l 



Pi c — a, 6; c — 1 



1 - z 
c[c- 1 + (a - 6)z] 



z - 1 
a(c — 6)z 
z-1 5 



2 Pi I c - a, 6; c; 



Pi c — a, 6; c + 1; 



z-1 



0. (C29) 



If we now c ombi ne this relationship with the linear trans- 
formation (C18), we obtain the following recurrence for- 
mula, where the first and the third parameter of the hy- 
pergeometric series change simultaneously: 

c(l - c) 2 F 1 (a- l,b;c— l;z) 
- c[c - 1 + (a — 6)z] 2 Pi(a, 6; c; z) 

+ a(c-6)z 2 Pi(a+l,6;c+l;z) = 0. (C30) 



If we now proceed in @, @, and in (jcToJ) - ( |cT5| ) in 

exactly the same way, we obtain the following recurrence 
formulas: 

c[c — a + (a — 6)z] 2 Pi (a, 6; c; z) 

— c(c — a) 2 Pi (a — 1, 6; c; z) 
-a(c-6)z 2 Pi(a+l,6;c+l;z) = 0, (C31) 

(1 - c) 2 Pi(a - 1, 6; c - 1; z) 
+ (a - l) 2 Pi(a, 6;c;z) 

+ (c-a) 2 Pi(a -1,6; c;z) = 0, (C32) 
c(l - z) 2 Pi(a + 1, 6; c; z) 

- c 2 Pi(a, 6;c;z) 

+ (c- 6)z 2 Pi(a+ l,6;c+ l;z) = 0, (C33) 
a(l — z) 2 Pi(a + 1, 6; c; z) 
+ (1 - c) 2 Pi(a- l,6;c- l;z) 

+ [c - a - 1 + (a - 6)z] 2 Pi(a, 6; c; z) = , (C34) 
c[b + (a — 6)z] 2Pi(o, 6; c; z) 

— 6c(l — z) 2 Pi(a, 6 + 1; c; z) 

- a(c- 6)z 2 Pi (a + 1,6; c+ l;z) = 0, (C35) 
(1 - c) 2 Pi(a- l,6;c- l;z) 
+ (c - 6 - 1) 2 Pi(a, 6; c; z) 

+ 6(l-z) 2 Pi(a,6+l;c;z) = 0, (C36) 
c 2 Pi(a, 6 - 1; c; z) 

- c 2 Pi(a, 6;c;z) 

+ az 2 F 1 (a + l,6;c+ l;z) = 0, (C37) 
(c - 6) 2 Pi(a, 6 - 1; c; z) 
+ (1 - c) 2 Pi(a- l,6;c- l;z) 

+ [6 - 1 + (a - 6)z] 2 Pi (a, 6; c; z) = . (C38) 

Next, we replace in (|C3| ) 6 by c — 6 and z by z/(z — 1). 
This yields: 



■ 2 Pi a,c-6;c- 1; 



1 - z 
c[c — 1 — (a — 6)z] 



z-1 
(c — a)6z 
z - 1 



z-1, 
2P1 I a, c - 6; c; 



z-1 



2^1 



Pi a, c — 6; c + 1: 



z-1 



0. (C39) 



If we now c ombi ne this relationship with the linear trans- 
formation (C17), we obtain the following recurrence for- 
mula, where the second and the third parameter of the 
hypergeometric series change simultaneously: 

c(l - c) 2 Pi(a, 6 - 1; c - 1; z) 

+ c[c - 1 - (a — 6)z] 2 Pi(a, 6; c; z) 

- (c-a)6z 2 Pi(a,6+l;c+l;z) = 0. (C40) 



If we now proceed in 0, @, and in ( |C10|) - (|GL5|) in 

exactly the same way, we obtain the following recurrence 
formulas: 



1G 



c[(a — b)z — a] 2-F1 (a, b; c; z) 
+ ac(l — z) 2-^1(0 + 1) b; c; z) 

+ (c — a)bz 2 F\(a, b + 1; c + 1; z) = 0, 
(1 - c) 2^1(0,6 -l;c-l;z) 
+ (c - a - 1) 2-F1 (a, 6; c; z) 

+ a(l — z) 2 F 1 (a + l,b;c; z) = 0, 
C2-Fi(o - 1, b; c; 2) 

— C2-Fi(a, 6; c; z) 

+ bz 2 F 1 (a, b + l;c + l;z) = 0, 
(c — a) 2-Fi(a — 1, b; c; z) 
+ (l-c) 2 Fi(a,&-l;c-l;;s) 

+ [a — 1 — (a — 6)z] 2-F1 (a, 6; c; z) = , 
c[c — b — (a — b)z] 2 F\{a, 6; c; z) 

— c(c- 6) 2-F1 (a, 6— l;c;z) 

— (c — a)bz 2-F1 (a, 6 + 1; c + 1; z) = 0, 
(l-c) 2 Fi(M-l;c-l;z) 
+ (b-l) 2 Fi(a,b;c;z) 

+ (c-b) 2 Fi(a,b-l;c;z) = 0, 
c(l - z) 2-Fi(a, 6 + 1; c; z) 

— C2-Fi(a, 6; c; z) 

+ (c-a)z 2-F1 (a, b+ l;c+ l;z) = 0, 
6(1 - 2)2-^1(0, 6 +l;c;z) 
+ (l-c) 2 Fi(a,b-l;c-l;z) 

+ [c — b — 1 — (a — b)z] 2-F1 (a, 6; c; z) = 



(C41) 



(C42) 



(C43) 



(C44) 



(C45) 



(C46) 



(C47) 



(C48) 



T he tw o groups of recursions flC30p - ( |C38|) and ( |C40| ) 
- ( |C48|) , respectively, are not independent. They can be 
transformed into each other by interchanging a and b. 
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TABLES 



TABLE I. Summation of the divergent hypergeometric series g2-Fi(l, 1; 2; —z) — ln(l + z) for z — 7/2 with the help of 



Wynn's epsilon algorithm, Brezinski's theta algorithm, and the Levin-type transformations (£,s n ) and <5j, n; (£, s n ) 
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TABLE II. Summation of the hypergeometric series 2 Fi(2/3, 4/3; 1/3; z) = (1 + z)(l - z)~ 5/3 for z = (1 + iv / 3)/2 with the 



help of Wynn's epsilon algorithm and Levin's transformation rf^,"'(C 







Sn 

Eq. ( 


(*) 
3.12 


) 










(n — 2[n 
-2[[n/2L 

Eq. (| 


/al) 
1) 












d< 0) (l, 
Eq. ( 


s (z)) 
A22|) 











1.00 


000 










1 


000 


000 


000 










1 


000 


000 


000 000 














1 


2.33 


333 


+ 


i2.30 


940 


2 


333 


333 


334 


+ 


i 2.309 


401 


077 


-0 


531 


737 


773 153 


+ 


il 


384 


198 


772 


334 


2 


0.38 


889 + 


i5.67 


728 


-0 


531 


737 


773 


+ 


i 1.384 


198 


772 


-1 


131 


921 


176 770 


+ 


il 


411 


935 


517 


843 


3 


-4.54 


938 


+ 


i5.67 


728 


-0 


953 


868 


768 


+ 


i 1.632 


539 


413 


-1 


120 


276 


444 364 


+ 


il 


322 


600 


394 


819 


4 


-7.49 


177 


+ 


i0.58 


091 


-1 


136 


352 


936 


+ 


i 1.326 


189 


971 


-1 


112 


622 


137 231 


+ 


il 


326 


484 


602 


288 


5 


-4.11 


180 




i5.27 


337 


-1 


115 


979 


051 


+ 


i 1.318 


232 


359 


-1 


113 


364 


776 960 


+ 


il 


326 


919 


577 


184 


6 


3.46 


967 




i5.27 


337 


-1 


111 


974 


302 


+ 


i 1.326 


730 


914 


-1 


113 


348 


742 255 


+ 


il 


326 


817 


605 


791 


7 


7.64 


993 


+ 


i 1.96 


705 


-1 


113 


209 


082 


+ 


i 1.327 


290 


958 


-1 


113 


338 


972 201 


+ 


il 


326 


827 


801 


234 


8 


3.09 


756 


+ 


i9.85 


197 


-1 


113 


429 


777 


+ 


i 1.326 


838 


273 


-1 


113 


340 


931 714 


+ 


il 


326 


828 


116 


649 


9 


-6.72 


205 


+ 


i9.85 


197 


-1 


113 


348 


679 


+ 


i 1.326 


799 


292 


-1 


113 


340 


813 708 


+ 


il 


326 


827 


863 


160 


10 


-11.9 


767 


+ 


i0.75 


063 


-1 


113 


334 


783 


+ 


i 1.326 


827 


026 


-1 


113 


340 


793 311 


+ 


il 


326 


827 


897 


184 


11 


-6.38 


818 




i8.92 


900 


-1 


113 


340 


291 


+ 


i 1.326 


829 


768 


-1 


113 


340 


798 965 


+ 


il 


326 


827 


896 


881 


12 


5.43 


724 




i8.92 


900 


-1 


113 


341 


213 


+ 


i 1.326 


827 


964 


-1 


113 


340 


798 476 


+ 


il 


326 


827 


896 


233 


13 


11.6 


655 


+ 


i 1.85 


859 


-1 


113 


340 


832 


+ 


i 1.326 


827 


770 


-1 


113 


340 


798 438 


+ 


il 


326 


827 


896 


343 


14 


5.12 


954 


+ 


i 13.1 


791 


-1 


113 


340 


767 


+ 


i 1.326 


827 


891 


-1 


113 


340 


798 455 


+ 


il 


326 


827 


896 


339 


15 


-8.54 


354 


+ 


il3.1 


791 


-1 


113 


340 


796 


+ 


i 1.326 


827 


905 


-1 


113 


340 


798 453 


+ 


il 


326 


827 


896 


338 


Exact 














-1.113 340 798 + 


i 1.326 827 896 


-1.113 340 798 453 


+ 


i 1.326 827 896 338 
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TABLE III. Acceleration of the convergence of the hypergeometric series 2-F\(3/7, 5/2; 7/2; z) for z = 77/100 with the help 
of Wynn's epsilon algorithm, Brezinski's theta algorithm, and the Levin-type transformations (£, s n ) and <k, (Ci s n)- 



«n(z) 



Eq. (4.3) 



(n-2[n/2]) 
£ 2[[n/2]l 



Eq. (2.5) 



q(n-3[[n/31) 
*2[[n/3|_ 



Eq. (A9) 



* 0) (l,so( Z )) 



Eq. (A.22) 



Eq. 



, s (z)) 

^1 






1 


000 


000 


000 


000 


1 


000 


000 


000 


000 


1 


000 


000 


000 


000 


1 


000 


000 


000 


000 


1 


000 


000 


000 


000 


1 

1 


1 


235 


1 14 


285 


71 A 

1 14 


1 


OO K 

235 


71 A 

1 14 


285 


71 A 

1 14 


1 


OO K 

235 


71 A 

1 14 


285 


71 A 

1 14 


1 


1 1 1 

411 


A07 
92 / 


77 

8/ / 


A A 7 
94/ 


1 


1 1 1 

411 


A07 

92/ 


77 

all 


A A 7 
94/ 


2 


1 


336 


547 


619 


048 


1 


411 


927 


877 


947 


1 


336 


547 


619 


048 


1 


448 


820 


807 


146 


1 


448 


820 


807 


146 


3 


1 


387 


972 


619 


048 


1 


441 


496 


598 


639 


1 


485 


041 


451 


400 


1 


459 


648 


832 


810 


1 


460 


767 


058 


064 


4 


1 


A~\ Q 
4l0 


cai 
091 


503 


003 


1 


A K7 

40 1 


A AA 

49U 


ceo 
562 


007 

32 1 


1 


A 7 1 

4/1 


Ton 
(89 


7P i 
/61 


A A 7 
44/ 


1 


1 c o 
4o2 


08/ 


OCA 

2bU 


952 


1 


4b3 


07 1 

2/1 


1 O K 

135 


352 


5 


1 


433 


666 


279 


063 


1 


461 


253 


523 


425 


1 


467 


349 


021 


834 


1 


463 


515 


006 


731 


1 


463 


725 


027 


496 


6 


1 


444 


100 


773 


353 


1 


463 


029 


920 


830 


1 


463 


655 


490 


309 


1 


463 


733 


041 


236 


1 


463 


795 


467 


606 


7 


1 


450 


702 


748 


500 


1 


463 


503 


143 


435 


1 


463 


758 


479 


471 


1 


463 


788 


732 


403 


1 


463 


805 


515 


965 


8 


1 


454 


973 


597 


660 


1 


463 


711 


233 


069 


1 


463 


790 


047 


799 


1 


463 


802 


617 


705 


1 


463 


806 


889 


463 


9 


1 


457 


785 


549 


506 


1 


463 


770 


364 


632 


1 


463 


806 


019 


078 


1 


463 


806 


018 


807 


1 


463 


807 


072 


295 


10 


1 


459 


663 


708 


382 


1 


463 


795 


262 


347 


1 


463 


806 


921 


084 


1 


463 


806 


841 


297 


1 


463 


807 


096 


230 


11 


1 


460 


933 


204 


660 


1 


463 


802 


626 


156 


1 


463 


807 


092 


160 


1 


463 


807 


038 


364 


1 


463 


807 


099 


326 


12 


1 


461 


799 


964 


188 


1 


463 


805 


637 


323 


1 


463 


807 


125 


890 


1 


463 


807 


085 


259 


1 


463 


807 


099 


724 


13 


1 


462 


396 


866 


847 


1 


463 


806 


552 


580 


1 


463 


807 


118 


461 


1 


463 


807 


096 


360 


1 


463 


807 


099 


774 


14 


1 


462 


811 


003 


597 


1 


463 


806 


919 


041 


1 


463 


807 


119 


752 


1 


463 


807 


098 


978 


1 


463 


807 


099 


781 


15 


1 


463 


100 


213 


038 


1 


463 


807 


032 


663 


1 


463 


807 


119 


214 


1 


463 


807 


099 


593 


1 


463 


807 


099 


781 


16 


1 


463 


303 


343 


588 


1 


463 


807 


077 


439 


1 


463 


807 


118 


928 


1 


463 


807 


099 


737 


1 


463 


807 


099 


782 



Exact 



1.463 807 099 782 



1.463 807 099 782 



1.463 807 099 782 



1.463 807 099 782 



1.463 807 099 782 
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TABLE IV. Acceleration of the convergence of the hypergeometric series 2-F\(3/7, 5/2; —7/2; 
of Aitken's iterated A 2 process, Wynn's epsilon algorithm, and Brezinski's theta algorithm. 



for 2 = 77/100 with the help 



»(*) 



A 



Eq. (4.8) 



(n-2[[n/2I) 



Eq. (A.4) 



(n-2[n/2j) 
e 2[[n/21 



Eq. (2.5) 



q(n-3In/31) 
7 2[[n/3L 



Eq. (A£) 






1 


000 


000 


10+° 


1 


000 


000 


000 


000 


10+° 


1 


000 


000 


000 


000 


10+° 


1.000 


000 


000 


000 


10+° 


1 


7 


642 


857 


lO" 1 


7 


642 


857 


142 


857 


10" 1 


7 


642 


857 


142 


857 


IO" 1 


7.642 


857 


142 


857 


10 


-1 


2 


9 


457 


857 


io- 1 


8 


668 


280 


871 


671 


nr 1 


8 


668 


280 


871 


671 


io- 1 


9.457 


857 


142 


857 


10 


-1 


3 


6 


063 


807 


10" 1 


8 


275 


261 


324 


042 


10" 1 


8 


275 


261 


324 


042 




8.485 


759 


131 


562 


10 


-1 


4 


3 


070 


461 


10+° 


8 


535 


660 


811 


732 


io- 1 


8 


609 


189 


507 


196 


nr 1 


8.494 


334 


311 


697 


10 


-1 


5 


2 


491 


700 


10+ 1 


8 


361 


725 


059 


741 


io- 1 


8 


157 


392 


107 


013 


nr 1 


8.453 


777 


105 


530 


10 


-1 


6 


1 


010 


158 


10+ 2 


8 


516 


151 


469 


394 


io- 1 


8 


896 


493 


143 


013 


nr 1 


8.487 


841 


750 


682 


10 


-1 


7 


2 


839 


789 


10+ 2 


8 


209 


992 


692 


837 


HP 1 


4 


152 


615 


385 


984 


io- 1 


8.486 


139 


776 


736 


10 


-1 


8 


6 


390 


582 


10+ 2 


8 


549 


782 


290 


518 


HT 1 


1 


979 


200 


475 


641 


10+° 


8.486 


753 


544 


589 


10 


-1 


9 


1 


236 


512 


10+ 3 


9 


170 


319 


462 


988 


HT 1 


1 


579 


725 


923 


915 


10+ 1 


8.486 


400 


009 


032 


10 


-1 


10 


2 


143 


447 


10+ 3 


8 


324 


631 


947 


362 


io- 1 


-9 


032 


991 


315 


812 


10+ 1 


8.486 


226 


769 


041 


10 


-1 


11 


3 


416 


644 


10+ 3 


7 


319 


750 


281 


906 


io- 1 


-4 


871 


173 


637 


209 


10+ 2 


8.486 


314 


552 


725 


10 


-1 


12 


5 


097 


263 


10+ 3 


8 


296 


736 


361 


657 


io- 1 


9 


550 


458 


201 


456 


10+ 3 


8.486 


265 


235 


534 


10 


-1 


13 


7 


207 


772 


10+ 3 


2 


786 


711 


790 


541 


10+° 


3 


152 


031 


922 


861 


10+ 4 


8.486 


238 


887 


649 


10 


-1 


14 


9 


751 


014 


10+ 3 


2 


417 


040 


069 


989 


10+° 


9 


017 


022 


343 


102 


10+ 4 


8.486 


245 


607 


105 


10 


-1 


15 


1 


271 


111 


10+ 4 


2 


660 


054 


094 


715 


10+° 


9 


807 


349 


664 


879 


10+ 4 


8.486 


241 


968 


013 


10 


-1 


16 


1 


605 


569 


1Q +4 


1 


583 


428 


788 


102 


10+° 


1 


006 


124 


741 


024 


10+ 5 


8.486 


240 


157 


205 


10 


-1 


17 


1 


973 


906 


10+ 4 


2 


773 


669 


575 


356 


10+° 


1 


009 


114 


496 


621 


10+ 5 


8.485 


741 


255 


417 


10 


-1 


IS 


2 


370 


573 


10+ 4 


2 


553 


560 


837 


044 


10+° 


1 


009 


964 


826 


770 


10+ 5 


8.486 


242 


524 


145 


10 


-1 


19 


2 


789 


405 


10+ 4 


2 


384 


039 


609 


516 


10+° 


1 


010 


111 


063 


638 


10+ 5 


8.486 


355 


858 


960 


10 


-1 


20 


3 


223 


965 


10+ 4 


2 


454 


155 


706 


579 


10+° 


1 


010 


153 


078 


626 


10+ 5 


8.492 


820 


131 


138 


10 


-1 


21 


3 


667 


836 


10 + 


2 


384 


021 


379 


340 


-i ,-,-1-0 

10 + 


1 


010 


161 


882 


819 


10 + 


8.486 


069 


295 


394 


10 


-1 


22 


4 


114 


878 


1Q +4 


2 


454 


149 


413 


545 


10+° 


1 


010 


164 


516 


908 


10+ 5 


8.470 


533 


544 


441 


10 


-1 


23 


4 


559 


414 


10+ 4 


7 


889 


100 


799 


808 


10+ 2 


1 


010 


165 


130 


151 


10+ 5 


2.015 


803 


318 


335 


10+° 


24 


4 


996 


370 


10+ 4 


7 


855 


297 


113 


412 


10+ 2 


1 


010 


165 


321 


636 


10+ 5 


8.519 


514 


381 


837 


10 


-1 


25 


5 


421 


360 


10+ 4 


7 


916 


396 


789 


254 


10+ 2 


1 


010 


165 


369 


332 


10+ 5 


1.185 


791 


071 


696 


10 


fO 


26 


5 


830 


726 


10+ 4 


7 


874 


382 


737 


052 


10+ 2 


1 


010 


165 


384 


712 


10+ 5 


-6.192 


361 


264 


079 


10 


fO 


27 


6 


221 


545 


10+ 4 


7 


780 


517 


151 


440 


10+ 2 


1 


010 


165 


388 


746 


10+ 5 


-2.153 


607 


969 


121 


10 


fO 


28 


6 


591 


599 


10+ 4 


7 


781 


235 


296 


899 


10+ 2 


1 


010 


165 


390 


075 


10+ 5 


-7.845 


642 


420 


775 


10+° 


29 


6 


939 


333 


10+ 4 


7 


780 


882 


446 


428 


10+ 2 


1 


010 


165 


390 


438 


10+ 5 


2.998 


151 


673 


644 


10+ 2 


30 


7 


263 


789 


10+ 4 


7 


781 


060 


228 


285 


10+ 2 


1 


010 


165 


390 


560 


10+ 5 


3.199 


681 


242 


234 


10 


H 



Exact 



1.010 165 



10 



+5 



1.010 165 390 611 • 10 



1.010 165 390 611 



10 



+5 



1.010 165 390 611 
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TABLE V. Acceleration of the convergence of the hypergeometric series 2^*1 (3/7, 5/2; —7/2; z) for z — 77/100 with the help 
of the iteration of Brezinski's theta algorithm and the Levin-type transformations (£, s n ) and 5t (£, s n ). 



s n (z) 



(n-3[[n/3]]) 



Eq. 



[n/31 



Eq. (A15) 



dk 0) (l,s (2)) 



Eq. 



5i 0) (Mo(*)) 
Eq. (S 






1.000 


000 


10+° 


1.000 


000 


000 


000 


10+° 


1.000 


000 


000 


000 


10+° 


1.000 


000 


000 


000 


10+° 


1 


7.642 


857 


10- 1 


7.642 


857 


142 


857 


nr 1 


8.668 


280 


871 


671 


10- 1 


8.668 


280 


871 


671 


10- 1 


2 


9.457 


857 


10- 1 


9.457 


857 


142 


857 


10- 1 


8.384 


394 


404 


076 


10- 1 


8.384 


394 


404 


076 


10- 1 


3 


6.063 


807 


nr 1 


8.485 


759 


131 


562 


10- 1 


8.564 


401 


754 


536 


10- 1 


8.543 


029 


016 


450 


10- 1 


4 


3.070 


461 


10+° 


8.494 


334 


311 


697 


10- 1 


8.410 


428 


387 


954 


10- 1 


8.451 


555 


013 


928 


10" 1 


5 


2.491 


700 


10+ 1 


8.453 


777 


105 


530 


10- 1 


8.581 


719 


847 


922 


10- 1 


8.511 


715 


495 


996 


10- 1 


6 


1.010 


158 


10+ 2 


8.487 


877 


036 


926 


10- 1 


8.339 


732 


531 


530 


10- 1 


8.466 


755 


128 


315 


HT 1 


7 


2.839 


789 


10+ 2 


8.485 


959 


518 


738 


10- 1 


8.765 


035 


290 


451 


10- 1 


8.504 


596 


561 


867 


10- 1 


8 


6.390 


582 


10+ 2 


8.489 


925 


522 


131 


10- 1 


7.853 


488 


289 


879 


w 1 


8.469 


220 


701 


568 


10- 1 


9 


1.236 


512 


10+ 3 


8.486 


664 


090 


804 


10- 1 


1.019 


118 


632 


891 


10+° 


8.505 


456 


725 


376 


10- 1 


10 


2.143 


447 


10+ 3 


8.335 


704 


458 


916 


10- 1 


3.141 


144 


923 


964 


10- 1 


8.465 


273 


890 


370 


lO" 1 


11 


3.416 


644 


10+ 3 


8.497 


498 


135 


494 


10- 1 


2.776 


012 


986 


506 


10+° 


8.513 


020 


606 


999 


10" 1 


12 


5.097 


263 


10+ 3 


8.501 


342 


364 


535 


10- 1 


-7.045 


936 


018 


531 


10+° 


8.452 


752 


921 


910 


10- 1 


13 


7.207 


772 


10+ 3 


8.498 


411 


613 


838 


10- 1 


3.717 


678 


472 


983 


10+ 1 


8.532 


975 


348 


760 


nr 1 


14 


9.751 


014 


10+ 3 


8.497 


632 


848 


220 


10- 1 


-1.858 


061 


858 


604 


10+ 2 


8.421 


055 


650 


264 


10- 1 


15 


1.271 


111 


10+ 4 


8.497 


961 


772 


597 


10- 1 


1.044 


967 


163 


943 


10+ 3 


8.583 


855 


953 


130 


lO" 1 


16 


1.605 


569 


10+ 4 


8.498 


474 


972 


833 


10- 1 


-7.005 


065 


553 


184 


10+ 3 


8.338 


036 


004 


332 


10- 1 


17 


1.973 


906 


10+ 4 


7.503 


910 


591 


903 


10- 1 


3.081 


110 


326 


042 


10+ 4 


8.721 


871 


069 


015 


lO" 1 


18 


2.370 


573 


10+ 4 


8.499 


855 


007 


426 


10- 1 


1.465 


687 


756 


463 


10+ 5 


8.104 


110 


955 


846 


10- 1 


19 


2.789 


405 


10+ 4 


7.414 


199 


029 


835 


10- 1 


9.719 


482 


708 


581 


10+ 4 


9.126 


025 


807 


834 


10- 1 


20 


3.223 


965 


1Q +4 


6.852 


823 


777 


347 


10- 1 


1.014 


869 


541 


810 


10+ 5 


7.392 


774 


724 


769 


10- 1 


21 


3.667 


836 


10+ 4 


6.800 


566 


895 


738 


10- 1 


1.009 


649 


257 


196 


10+ 5 


1.040 


040 


245 


220 


10+° 


22 


4.114 


878 


10+ 4 


6.778 


745 


199 


478 


10" 1 


1.010 


218 


868 


540 


10+ 5 


5.070 


943 


906 


098 


KT 1 


23 


4.559 


414 


10+ 4 


6.776 


231 


440 


878 


10- 1 


1.010 


160 


173 


901 


10+ 5 


1.469 


832 


601 


187 


10+° 


24 


4.996 


370 


10+ 4 


6.777 


792 


652 


237 


10- 1 


1.010 


165 


871 


479 


10+ 5 


-3.005 


355 


825 


647 


10- 1 


25 


5.421 


360 


10+ 4 


6.807 


217 


929 


960 


10- 1 


1.010 


165 


348 


613 


10+ 5 


3.008 


947 


036 


607 


10+° 


26 


5.830 


726 


1Q +4 


6.775 


201 


071 


781 


10- 1 


1.010 


165 


394 


094 


10+ 5 


-3.273 


286 


658 


302 


10+° 


27 


6.221 


545 


10+ 4 


6.796 


867 


331 


534 


10- 1 


1.010 


165 


390 


335 


10+ 5 


8.821 


765 


469 


703 


10+° 


28 


6.591 


599 


10+ 4 


1.050 


749 


944 


450 


10+° 


1.010 


165 


390 


631 


10+ 5 


-1.477 


966 


683 


169 


10+ 1 


29 


6.939 


333 


10+ 4 


-6.085 


082 


692 


546 


10- 1 


1.010 


165 


390 


609 


10+ 5 


3.183 


343 


913 


267 


10+ 1 


30 


7.263 


789 


10+ 4 


-6.037 


663 


167 


562 


10- 1 


1.010 


165 


390 


611 


10+ B 


-6.136 


838 


085 


955 


10+ 1 


Exact 


1.010 165 


10+ 6 


1.010 165 390 611 


10+ 5 


1.010 165 390 611 


10+ 5 


1.010 165 390 611 


10+ 5 
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TABLE VI. Acceleration of the convergence of the hypergeometric series 2-Fi(3/7, 5/2; —7/2; z) for z = 77/100 with the help 
of Wynn's epsilon algorithm and the Levin-type transformations d^(^,s n ) and 8^(1^, s n ). 



D (22) 



Eq. (1.11) 



(„-2[n/2]) 
e 2[[n/2]l 



Eq. (2.5) 



4 0) (i, 

Eq. 



4 22) (-)) 



0A2 



^ 0) (1, 
Eq. 



•? 2) M) 



(A25) 






4 


114 


878 


10 4 


4 


114 


877 


834 


479 


10 4 


4 


114 


877 


834 


479 


10 4 
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1 
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14 


8 
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017 
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010 
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8 
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9 
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1 


010 


165 


390 


507 
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1 
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165 


339 
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1 


010 


165 


390 
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9 


302 


057 
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010 


165 


390 
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10 5 


1 


010 


165 


388 


497 


10 5 


1 


010 


165 


390 


609 


10 5 


19 


9 


408 


137 


10 4 


1 


010 


165 


390 


607 


10 5 


1 


010 


165 


393 


632 


10 5 


1 


010 


165 


390 


611 


10 5 


20 


9 


501 


598 


10 4 


1 


010 


165 


390 


610 


10 5 


1 


010 


165 


392 


321 


10 5 


1 


010 


165 


390 


611 


10 5 



Exact 



1.010 165 



10 5 



1.010 165 390 611 



10 5 



1.010 165 390 611 



10" 



1.010 165 390 611 



23 



(3) 

TABLE VII. Selected coefficients c„ of the renormalized 



ratios Ci 3) defined by Eq. (|B12|) 



weak coupling expansion (B6) as well as the corresponding 



n 




c (3) 




M3) 





1.000 


000 000 




-0.49218 


1 


-3.333 


333 333 


lO" 1 


-1.51799 


2 


-9.074 


074 074 


10~ 2 


0.32773 


3 


3.451 


646 091 


HT 1 


0.34954 


4 


-3.064 


808 585 




0.44181 


5 


4.145 


321 167 


10+ 1 


0.51354 


6 


-8.011 


680 849 


10+ 2 


0.57041 


7 


2.103 


995 759 


10+ 04 


0.61602 


8 


-7.225 


346 394 


10+ 05 


0.65320 


9 


3.148 


105 306 


10+ 07 


0.68398 


10 


-1.698 


432 299 


10+ 09 


0.70986 


11 


1.112 


192 278 


10+ 11 


0.73188 


12 


-8.693 


791 326 


10+ 12 


0.75085 


13 


7.998 


709 458 


10+ 14 


0.76734 


14 


-8.558 


133 512 


10+ 16 


0.78181 


15 


1.053 


809 185 


10+ 19 


0.79460 


25 


3.966 


243 637 


10+ 42 


0.87072 


50 


—8.551 


437 639 


10 +114 


0.93301 


75 


4.053 


097 428 


10 +198 


0.95481 


100 


-4.540 


614 326 


10 +289 


0.96591 


125 


1.552 


803 192 


10+ 386 


0.97263 


150 


-1.248 


762 212 


10 +487 


0.97714 


175 


4.338 


412 567 


10+ 591 


0.98037 


200 


-1.952 


374 463 


ig+699 


0.98280 


225 


4.623 


901 073 


10+ 809 


0.98470 


250 


-2.864 


322 271 


10 +922 


0.98622 


275 


2.654 


720 800 


10 +1037 


0.98746 


300 


-2.331 


943 009 


10 +1154 


0.98850 
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TABLE VIII. Infinite coupling limit of the sextic an- 



harmonic oscillator. 



Convergence of the approximants n, 3 



(n,i) 



defined by Eq. (B1J) with the three highest possible values of 
n < 299 - I for < I < 12. 



/ 


n 







297 


1.144 802 453 797 051 992 831 




298 


1.144 802 453 797 052 042 002 




299 


1.144 802 453 797 052 088 337 


1 


296 


1.144 802 453 797 052 586 551 




297 


1.144 802 453 797 052 601 411 




298 


1.144 802 453 797 052 615 225 


2 


295 


1.144 802 453 797 052 759 384 




296 


1.144 802 453 797 052 761 860 




297 


1.144 802 453 797 052 764 077 


3 


294 


1.144 802 453 797 052 783 941 




295 


1.144 802 453 797 052 783 501 




296 


1.144 802 453 797 052 783 043 


4 


293 


1.144 802 453 797 052 775 728 




294 


1.144 802 453 797 052 775 182 




295 


1.144 802 453 797 052 774 657 


5 


292 


1.144 802 453 797 052 767 958 




293 


1.144 802 453 797 052 767 714 




294 


1.144 802 453 797 052 767 482 


6 


291 


1.144 802 453 797 052 764 644 




292 


1.144 802 453 797 052 764 570 




293 


1.144 802 453 797 052 764 501 


7 


290 


1.144 802 453 797 052 763 692 




291 


1.144 802 453 797 052 763 725 




292 


1.144 802 453 797 052 763 758 


8 


289 


1.144 802 453 797 052 765 158 




290 


1.144 802 453 797 052 764 915 




291 


1.144 802 453 797 052 764 334 


9 


288 


1.144 802 453 797 052 746 941 




289 


1.144 802 453 797 052 719 382 




290 


1.144 802 453 797 052 727 384 


10 


287 


1.144 802 453 797 049 918 777 




288 


1.144 802 453 797 053 441 988 




289 


1.144 802 453 797 055 805 549 


11 


286 


1.144 802 453 797 422 201 710 




287 


1.144 802 453 797 325 472 977 




288 


1.144 802 453 796 999 625 608 


12 


285 


1.144 802 453 789 565 166 506 




286 


1.144 802 453 760 484 424 039 




287 


1.144 802 453 769 883 065 332 
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